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PREFACE 


Our  primary  objective  is  to  determine  the  space  charge  dis¬ 
tribution  surrounding  a  high-speed  vehicle  at  high  altitudes  fol¬ 
lowing  an  atomic  bursts  To  date,  this  problem  remains  unsolved. 

Indeed,  no  procedure,  theoretical  or  experimental,  is  available 
for  the  simpler  problem  of  flight  not  involving  an  atomic  burst 
Although  references  1,  2,  and  3  describe  an  instrument  to  measure 
the  undisturbed  free  stream  ionization  for  argon  gas  in  a  low- 
density  supersonic  plasma  jet  stream,  the  assumptions  of  this  theory 
do  not  appear  valid  for  our  use  in  air. 

To  begin  a  solution  of  the  above  problem,  we  decided  to  analyze 
in  exact  detail  a  much  simpler  problem — one  from  which  we  could  be¬ 
gin  to  build  a  theory  and  readily  confirm  experimentally  as  we  pro¬ 
ceeded.  Such  a  problem  is  presented  in  this  report.  Here,  we  show 
a  theoretical  procedure  for  solving  the  ion-electron  distribution 
in  the  stagnation- point  boundary  layer  of  a  blunt  body  and  confirm 
this  procedure  with  measurements  in  a  shock  tube.  This  problem  is 
solved  exactly  with  the  use  of  electronic  computers  by  simultaneously 
solving  the  gasdynamic  and  electric  field  equations. 

The  neutral  gas  equations  for  a  dissociated  gas  cannot  be  solved 
in  a  general  way.  At  best,  solutions  can  be  obtained  for  the  stagna¬ 
tion  region  of  a  blunt  body  behind  the  detached  shcck  and  for  conical 
and  flat  plate  bodies.  This  limits  the  type  of  body  shapes  for  which 
we  can  expect  to  obtain  solutions.  In  addition,  we  will  be  limited 
by  the  precision  to  which  we  can  describe  the  ionization  and  dissocia¬ 
tion  effects  on  the  neutral  gas  equations,  and  these  will  be  dependent 
on  the  special  conditions  for  the  case  considered. 

The  solution  of  the  blunt  body  stagnation-point  problem  is  of 
special  interest  in  its  application  as  a  device  for  the  determination 
of  the  space  charge  distribution  about  bodies  of  an  arbitrary  shape. 
That  is,  we  propose  to  mount  such  a  probe  or  probes  on  a  vehicle.  The 
probe  measurements  will  then  be  interpreted  in  terms  of  the  space  charge 
at  a  point  in  a  region  undisturbed  by  the  probe.  In  this  way,  the 
blunt  body  stagnation-point  probes  can  be  used  as  instruments  to  map 
the  space  charge  distribution  about  a  vehicle  of  arbitrary  shape.  In 
turn,  this  procedure  may  have  direct  application  to  our  primary  prob¬ 
lem  involving  an  atomic  burst. 

To  investigate  the  above  possibilities  further,  our  immediate 
plan  includes  the  extension  of  the  theory  and  measurements  for  the 
blunt  body  stagnation-point  region  given  in  the  present  paper  to  higher 
flow  Mach  numbers,  gas  velocities,  and  altitudes.  For  these  experi¬ 
mental  measurements,  we  plan  the  use  of  a  shock  tunnel.  Then,  depend¬ 
ing  on  the  success  of  this  preliminary  shock  tunnel  work,  we  will  con¬ 
sider  the  insertion  of  simple  body  shapes  in  the  shock  tunnel  and 


attempt  to  interpret  our  measurements  on  blunt  body  probes  in  terns 
of  the  ion-electron  density  at  a  point  im  the  neighborhood  of  the 
body.  Then,  by  a  point -by -point  determination,  we  hope  to  map  the 
charge  distribution  about  a  vehicle. 
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ABSTRACT 


A  model  of  the  ion-electron  flow  in  the  stagnation  region 
between  the  detached  shock  and  the  wall  of  a  Langmuir-type  probe 
is  proposed  for  ionized  air  in  a  shock  tube  at  ion  mass  fractions 
less  than  IQ"4.  Based  on  the  close  agreement  between  numerical 
calculations  and  measurements  of  the  probe  well  ion  current  den¬ 
sity,  the  proposed  model  gives  an  accurate  description  of  the  ion- 
electron  flow  in  the  stagnation-point  boundary  layer  when  negative 
potentials  are  applied  at  the  stagnation  point  of  the  probe.  The 
procedure  is  indicated  to  be  applicable  also  in  the  region  between 
the  detached  probe  shock  and  the  edge  of  the  boundary  layer.  The 
Ion  and  electron  densities  in  the  boundary  layer  are  shown  to  be 
strongly  dependent  on  the  equilibrium  conditions  for  the  high- tem¬ 
perature  gas  behind  the  shock.  Hence,  the  ion  and  electron  number 
densities  in  the  undisturbed  incident  flow  ahead  of  the  detached 
shock  cannot  be  obtained  unless  it  is  possible  to  calculate  the 
changes  in  the  number  densities  across  the  shock. 

The  measurements  were  made  using  bottled  dry  air  as  the  test 
gas  for  shock  tube  Mach  numbers  from  7.2  to  11.0,  electrode  poten¬ 
tials  (referred  to  the  plasms  potential)  of  -6.1,  -11.5  and  -22.3  v, 
initial  shock  tube  pressure  of  1  cm  Hg,  and  initial  shock  tube  tem¬ 
peratures  between  295  and  300°K.  The  ion  number  density  distribu¬ 
tion  through  the  boundary  layer  including  the  sheath  is  calculated 
for  an  electrode  potential  of  -11.5  v  at  shock  tube  Mach  numbers  of 
7,43,  9,88,  and  10.43.  Additional  calculations  are  presented  for 
the  determination  of  the  electric  field,  wall  ion  current  density, 
and  sheath  thickness  for  potentials  in  the  range  -2  to  -24  v.  The 
effects  of  some  of  the  lesser-known  parameters  on  the  ion  current 
density  calculations  are  also  indicated.  These  parameters  include 
the  Lewis  number  for  the  nitric  oxide  ion-electron  pair  in  the  par¬ 
tially  ionized  boundary  layer  flow,  initial  ion  number  density  at 
the  edge  of  the  boundary  layer,  stagnation-point  velocity  gradient, 
and  the  ion-electron  formation- recombination  rates  in  the  boundary 
layer. 
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(J>  sheath  potential  referred,  to  plasma,  v 
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1. 


INTRODUCTION 


It  was  recently  shown  how  probes  could  be  used  in  a  low  den- 


Stagnation-point  probe 


incident  ion  and  electron  densities  (ref  1,  2,  3,).  As  shown  in 
the  above  sketch,  the  probe  consists  of  a  flat-headed  circular 
cylinder  with  an  electrode  placed  at  the  stagnation  point  of  the 
incident  flow.  For  air  at  equilibrium  fractional  ionization  densi¬ 
ties  less  than  10“4,  all  electrons  and  ions  may  be  assumed  formed 
through  the  reaction  (ref  4) 

N  +  0  £  N0+  +  e  (1.1) 


To  measure  ion  current  only,  the  electrode  is  biased  negatively  with 
respect  to  the  plasma.  The  negative  potential  causes  an  ion  sheath 
to  form  at  (the  electrode  and,  if  sufficiently  strong,  will  repel, 
all  electrons.  A  boundary  layer  is  formed  at  the  probe  wall  and  its 
thickness  in  the  stagnation  region  is  much  less  than  the  probe  shock 
detachment  distance.  The  sheath  is  much  thinner  than  the  boundary 
layer.  Consequently,  outside  the  sheath,  ions  and  electrons  in  the 
boundary  layer  diffuse  through  the  gas  mixture  at  a  speed  corres¬ 
ponding  to  ambipolar  diffusion.  In  the  sheath  there  occurs  diffu¬ 
sion  of  ions. 

The  well-known  blunt-body  stagnation  point  solutions  for  the 
continuum  flow  of  a  nonionized  gas  are  applicable  for  slightly  ion¬ 
ized  gases.  These  solutions  give  the  enthalpy,  pressure,  -and  veloc¬ 
ity  distributions  for  the  gas  mixture  in  the  stagnation  region 
between  the  detached  shock  and  the  probe  wall.  The  equation  for  the 
electric  potential  together  with  the  conservation  equations  for  ions 
and  electrons  yields  solutions  for  the  electric  potential,  ion  and 
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electron  distributions  in  the  boundary  layer  between  the  shock  and 
the  probe  wall.  The  flow  and  fluid  properties  of  the  gas, mixture 
as  well  as  formation  and  recombination  rate  constants  of  the  various 
sir  constituents  are  known  in  the  entire  region  behind  the  shock. 

/  cordingly,  ion  and  electron  distributions  between  the  probe  wall 
and  detached  shock  can  bo  determined.  Moreover,  the  ion-electron 
density  of  the  incident  flow  may  be  obtained  when  the  variation  of 
this  quantity  across  the  shock  is  known.  For  neutral  plasma  flow 
behind  the  probe  shock,  due  to  ambipolar  diffusion,  the  electron  and 
ion  flows  and  density  clist  ri  but  ions  are  the  same,  excluding  the  thin 
sheath  formed  at  the  stagnation  point. 

For  the  low  gas  densities  in  the  plasma  jet,  frozen  flow  is 
assumed  in  the  entire  region  between  the  detached  probe  shock  and 
the  stagnation  point.  The  ion  sheath  is  a  small  fraction  of  the 
boundary  layer  thickness,  bince  the  fractional  ionization  is  less 
than  10“3  and  the  sheath  is  so  thin,  the  charged  particles  are  as¬ 
sumed  to  have  no  effect  on  the  enthalpy  and  motion  of  the  neutral  gas 
The  mean  free  path  of  the  gas  is  sufficiently  small  to  permit  the  use 
of  the  hypersonic  boundary  layer  equations.  Therefore,  outside  the 
sheath,  the  ion  motion  is  assumed  to  be  the  sum  of  the  gas  motion  and 
the  relative  motion  due  to  ion-electron  ambipolar  diffusion. 

In  a  low-density  plasma  jet,  the  ion  mean  free  path  is  larger 
than  the  sheath  thickness,  and  consequently,  ion  collisions  do  not 
occur  within  the  sheath.  Moreover,  since  the  ion  drift  velocity  is 
small  compared  with  the  ion  thermal  velocity,  the  potential  distribu¬ 
tion  within  the  sheath  is  given  by  the  familiar  Child’s  law, 


The  current  density  is  assumed  constant  within  the  sheath  by  virtue 
of  the  assumed  free-fall  flow  of  the  ions  through  the  sheath.  More¬ 
over,  J  is  independent  of  (j>w  and  within  the  sheath  (j>  oc  yV5  .  Thus, 
when  the  ion  density  ratio  across  the  probe  shock  is  known,  which 
here  is  assumed  equal  to  the  gas  density  ratio,  the  undisturbed  free 
stream  ion  or  electron  number  density  at  a  point  ahead  of  the  probe 
may  be  determined  from  the  measured  probe  current  density  by  solving 
the  hypersonic  boundary  layer  equations  and  the  sheath  equation. 

The  present  effort  is  an  extension  of  the  above  work  in  that 
ion-electron  densities  in  a  neutral  plasma  are  measured  in  a  shock 
tube  for  which  the  gas  is  about  103  times  more  dense.  This  increased 
density  changes  the  physics  of  the  flow  in  a  number  of  ways .  The 
region  between  the  detached  probe  shock  and  the  stagnation  point  is 
no  longer  frozen.  The  ion  number  density  increases  by  a  factor  of 
more  than  100  between  the  shock  and  the  edge  of  the  thermal  boundary 
layer.  As  shown  by  the  calculations  in  section  3.4,  the  ionization 
attains  thermodynamic  equilibrium  in  this  region. 
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Ion-electron  recombination  occurs  within  the  boundary  layer  up  to 
the  edge  of  the  sheath,  where  electrons  are  sharply  repelled  because 
of  their  low  energy  (approximately  0.3  v) .  For  a  probe  in  a  shock 
tube,  the  sheath  is  also  much  thinner  than  the  boundary  layer  and  the 
fractional  ionization  encountered  in  the  present  test  measurements 
and  calculations  was  always  less  than  10”4.  However,  the  flow  inci¬ 
dent  to  the  probe  is  now  supersonic,  and  the  corresponding  boundary 
layer  equations  are  applicable.  Outside  the  sheath,  the  ion  motion  is 
again  the  sum  of  the  gas  motion  and  the  relative  motion  due  to  ion- 
electron  ambipolar  diffusion.  Here,  the  diffusion  velocity  is  appre¬ 
ciably  affected  by  formation  and  recombination. 

Additional  ion-electron  formation  within  the  sheath  is  negligi¬ 
ble  when  the  sheath  is  a  small  fraction  of  the  boundary  layer.  Ac¬ 
cordingly,  formation  within  the  sheath  becomes  important  at  low  shock 
tube  Mach  numbers  and  large  probe  potentials. 

Inside  the  sheath,  ion  collisions  occur  mainly  with  neutral  parti¬ 
cles.  Hence,  the  ion  flow  through  the  sheath  is  largely  governed  by 
the  mobility  and  the  electric  field  resulting  from  the  potential  ap¬ 
plied  at  the  probe.  In  addition,  the  sheath  equation  includes  terms 
to  account  for  the  initial  ion  drift  velocity  across  the  sheath  due 
to  the  neutral  gas  motion  and  diffusion  caused  by  a  concentration 
gradient.  Here,  ion  thermal  diffusion  does  not  affect  the  probe  cur¬ 
rent  density.  The  current  density  increases  slightly  with  distance 
from  the  stagnation  point,  because  ions  are  carried  away  by  the  neutral 
gas  flow  tangential  to  the  probe  surface. 

Thus,  in  the  shock  tube,  a  determination  of  the  ion  or  electron 
number  density  distribution  in  the  stagnation  region  of  a  blunt  body 
behind  the  detached  shock  requires  the  solution  of  the  supersonic 
boundary  layer  equations  and  also  a  different  sheath  equation.  In  ad¬ 
dition,  account  must  be  taken  of  the  ion-electron  formation  and  recom¬ 
bination  processes,  which  in  turn  require  a  knowledge  of  the  local 
thermodynamic  conditions  and  the  rate  constants.  To  obtain  the  undis¬ 
turbed  incident  ion  or  electron  number  density,  the  variation  in  the 
number  density  across  the  probe  shock  must  also  be  known.  Since  the 
flow  is  not  frozen,  this  variation  is  not  determined  by  the  gas  den¬ 
sity  ratio  across  the  shock,  as  was  the  case  for  the  plasma  jet. 

The  following  presentation  describes  the  ion-electron  flow  in 
the  boundary  layer  behind  the  detached  probe  shock  for  the  shock  tube 
Mach  numbers  7.43,  9  88,  and  10.43  at  an  initial  shock  tube  pressure 
of  1  cm  Hg  and  temperature  of  300°  K  for  probe  potentials  in  the  range 
-5.4  to  -21.6  v  with  air  as  test  gas.  A  description  of  the  test  equip¬ 
ment  and  probe  measurements  of  the  stagnation-point  ion  current  density 
made  in  the  HDL  2  in.  by  2  in.  shock  tube  is  given  in  section  2. 

A  model  describing  the  physical  processes  involved  in  the  flow 
over  the  probes  permits  the  assumption  of  distinct  flow  regimes,  see 
section  3.  Using  this  model,  ion-electron  number  density  distributions 
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in  the  boundary  layer  behind  the  probe  shock  are  calculated  in  sec¬ 
tion  4.  The  validity  of  the  calculations  is  indicated  by  the  com¬ 
parison  between  the  calculated  and  experimental  values  for  the 
stagnation-point  ion  current  density. 

Calculations  are  included  in  section  4  to  show  the  effects  on 
the  ion  current  density  of  variation  of  some  of  the  lesser  known 
parameters.  The  parameters  varied  include  the  ratio  of  the  Prandtl 
to  Lewis  numbers  (the  Schmidt  number)  for  the  ion-electron  pairs  in 
the  partially  ionized  boundary  layer  flow,  initial  ion  number  den¬ 
sity  at  the  edge  of  the  boundary  layer,  stagnation  point  velocity 
gradient,  and  the  ion-electron  formation  and  recombination  rates  in 
the  boundary  layer. 

In  appendix  A,  for  the  described  test  conditions.  Coulomb  forces 
are  shown  to  have  a  negligible  effect  on  the  ion  mean  free  path  and 
the  coefficients  for  diffusion  and  viscosity. 

2 ♦  EXPERIMENTAL  MEASUREMENTS 

2 . 1  Test  Equipment 


The  test  measurements  were  made  in  the  HDL  2  in.  by  2  in. 
shock  tube  facility  The  unit  is  an  ordinary  single-diaphragm  straight 
shock  tube  in  which  the  shock  velocity  is  varied  by  controlling  the 
driver  gas  enthalpy;  the  shock  tube  Mach  number  was  varied  through  con¬ 
trolled  heat  audition  to  helium  driver  gas.  The  shock  tube  Mach  num¬ 
ber  is  determined  with  a  precision  of  better  than  1  percent,  using  six 
pressure  pickups  at  1-ft  intervals  along  the  tube.  All  the  reported 
data  were  obtained  using  bottled  dry  air  as  the  test  gas  at  an  initial 
pressure  of  1  cm  Hg  and  initial  temperatures  between  295  and  300°  K. 

The  shock  tube  comprises  two  sections,  one  of  which  is  mov¬ 
able  and  permits  the  lengthening  of  the  tube  for  inserting  test  sec¬ 
tions  of  lengths  up  to  30  in .  A  special  test  section  to  mount  probes 
was  constructed  from  a  single  piece  of  2- in.  thick  plexiglas  and  lo¬ 
cated  17  ft  from  the  aluminum  diaphragm.  Figures  1,  2,  and  3  show  the 
position  of  the  test  section  in  the  shock  tube  and  the  installation  of 
the  probes  with  axes  parallel  to  the  tube  centerline  In  figures  1, 

2,  and  4,  a  second  plexiglas  section  is  shown  adjacent  to  the  test  sec¬ 
tion  to  help  illustrate  the  probe  installation  and  flow  of  the  hot  gas 
over  the  probes  This  section  was  removed  during  the  test  measure¬ 
ments  to  obtain  a  flush  alignment  between  the  walls  of  the  test  sec¬ 
tion  and  the  shock  tube.  The  stagnation  point  of  all  the  probes  tested 
was  always  located  upstream  of  the  test  section-shock  tube  joint. 

Four  base  supports  were  placed  symmetrically  in  the  test 
section  with  the  centerline  of  each  support  parallel  to  and  1/2  in. 
from  the  vertical  centerline  of  the  shock  tube  (fig.  2  and  3)  This 
allowed  the  stagnation  point  of  each  of  four  probes  to  be  located  1/2 
i  ,  from  the  r.wo  near  walls  of  the  shock  uuhe  in  the  same  plane  normal 


to  the  shock  tube  axis.  Each  support  has  a  metal  insert  connected 
to  a  center  cable  and  to  the  center  electrode  of  an  external  plug, 
all  of  which  is  electrically  insulated  from  the  exposed  outer  surface 
of  the  base  support.  The  latter  was  electrically  connected  to  ground 
through  the  external  plug,  using  the  shock  tube  wall  as  ground.  A 
probe  was  threaded  into  the  metal  insert  of  each  support,  thus  con¬ 
necting  the  probe  electrode  to  the  external  plug.  A  teflon  washer 
was  used  to  insure  electrical  insulation  of  the  probe  electrode  from 
the  plasma  and  base  support  (fig  2)  The  probe  electrode  was  in¬ 
sulated  everywhere  from  the  plasma,  except  at  the  stagnation  point. 

Four  simultaneous  probe  current  measurements  were  recorded 
using  two  Tektronix  555  dual-beam  oscilloscopes.  Figure  4  shows  the 
flow  over  the  probes  during  a  typical  test,  the  radiation  of  the  hot 
gas  causing  the  film  exposure. 

A  total  of  five  sets  of  probes  was  constructed.  Three  sets 
were  flat-headed  circular  cylinders  each  of  2-in.  length  (excluding 
thread)  having  glass  insulation  over  the  electrode  with  the  following 
dimensions : 

(1)  0.02-in.  electrode,  0.25-in.  diameter  cylinder, 

(2)  0.06-in  jeledtrode,  0I25-in  diamdter  pylirider, 

(3)  0.02~in.  electrode,  0.10-in,  diameter  cylinder. 

The  third  set  was  broadened  to  0.25  in  diameter  at  the  base  in  order 
to  secure  firmly  in  the  support  (fig.  3) 

A  fourth  set  of  1-in.  length  probes  with  dimensions  otherwise 
identical  with  set  2  was  constructed  with  a  brass  sleeve  concentric 
with  the  cylinder  axis  and  running  the  entire  length  of  the  probe.  An 
epoxy  resin  insulator  was  used  between  the  sleeve  and  the  brass  elec¬ 
trode.  Two  of  these  probes  are  shown  in  figure  3,  outside  the  test 
section . 


In  the  fifth  set,  each  probe  consisted  of  a  conical  nose  with 
a  20-deg  vertex  angle  and  0.25-in  cylindrical  afterbody.  Excluding 
thread,  the  probes  were  of  2-in„  lengths.  The  electrode  was  exposed 
to  the  plasma  at  the  nose  tip  for  0.142  in.,  which  corresponds  to  a 
diameter  of  0.050  in.  Glass  insulation  covered  the  remainder  of  the 
probb . 


For  reasons  given  in  section  2.2,  current  measurements  are 
reported  only  for  the  first  three  sets  of  probes.  As  discussed  in 
section  2.2,  the  plasma  potential  may  be  assumed  to  be  about  +0.7  v, 
referred  to  the  shock  tube  wall.  Ion  currents  were  measured  by  using 
probe  potentials  of  -5.4,  -10.8,  and  -21.6  v,  with  the  shock  tube 
wall  as  ground.  These  potentials  were  obtained  by  placing  mercury 
cells  in  series,  each  cell  having  a  potential  of  1  35  v  and  capable  of 
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delivering  10  raa.  The  current  was  always  observed  to  be  less  than 
10  ma  and  the  steady'  state  value  never  exceeded  3  ma.  A  resistor 
was  placed  in  series  between  the  battery  and  ground,  and  the  oscillo¬ 
scope  was  connected  across  the  resistor.  The  resistor  had  a  value  of 
100  or  1000  ohms  and  was  always  used  so  that  the  steady-state  poten- 
tial  drop  across  it  and  thereby  the  reduction  in  the  probe  potential 
was  less  than  2  percent  of  the  applied  potential.  The  complete  cir¬ 
cuit  from  probe  to  ground  had  a  capacitance  of  4  x  10"10  f,  so  that 
the  time  constant  with  the  1000-ohm  resistor  was  0.4  x  10"6  sec"’1  . 

The  d-c  conductivity  of  the  undisturbed  air  stream  in¬ 
creases  with  increasing  Ms  and  amounts  to  about  0.01  mho/cm  at  Mg  = 
10.5  (fig.  5  of  ref  6).  Thus,  for  sufficiently  large  Ms,  it  was 
conceivable  that  the  gas  could  provide  a  conducting  path  between  ad¬ 
jacent  probes  or,  because  of  the  increased  conductivity  behind  the 
detached  probe  shock,  a  conducting  path  could  be  established  between 
a  probe  and  the  shock  tube  wall.  To  determine  whether  this  effect 
was  present,  at  the  higher  Mg,  current  density  measurements  were  ob¬ 
tained  using  both  100-  and  1000-ohm  resistors.  In  addition,  measure¬ 
ments  obtained  by  setting  all  four  probes  at  the  same  potential  were 
compared  at  the  same  potential  and  Mg  with  measurements  obtained  when 
each  of  the  four  probes  was  set  at  four  different  potentials .  In 
each  instance,  no  observable  effect  was  noted,  and  thus  it  appears 
that  the  gas  conductivity  did  not  influence  the  measurements. 

2.2  Test  Results 


Initially,  measurements  were  made  using  four  identical  probes 
biased  at  the  same  potential  to  observe  the  uniformity  of  the  ionized 
flow  across  the  tube  and  determine  bias  and  repeatability  of  the  probes 
As  shown  in  figure  5  for  (j)p  =  -10.8  v  on  the  0.06-in.  electrode,  0.25- 
in.  diameter  probes,  the  total  variation  of  the  four  measured  probe  cur 
rents  on  a  given  test  is  within  20  percent  and  is  usually  closer  to  10 
percent.  This  variation  does  not  appear  to  be  due  to  individual  pecu¬ 
liarities  in  the  probes,  although  the  results  of  the  measurements  were 
altered  when  the  probe  faces  were  not  clean  and  polished,  and  when 
occasionally  the  probes  were  deformed  during  a  test  series  by  diaphragm 
debris  (in  the  form  of  very  fine  aluminum  dust  particles  impacting  the 
probe  face).  After  about  20  firings,  the  currents  measured  by  one  or 
two  probes  would  suddenly  increase  50  per  cent  or  more  above  the  cur¬ 
rents  measured  by  the  remaining  probes.  This  was  apparently  due  to 
diaphragm  dust  that  had  coated  the  probe  face,  for  on  polishing  the 
probe  surface  with  a  fine  abrasive  the  probe  current  measurements  were 
restored  to  their  previous  values.  This  dust  coating  was  observable 
with  the  naked  eye  and  serveral  of  the  larger  flakes  along  the  test 
section  wall  were  noted,  by  means  of  an  ohmmeter,  to  be  an  electrical 
conductor. 

As  may  be  seen  from  figure  5,  an  error  of  1  percent  in  the 
determination  of  Ms  will  result  in  a  10-percent  error  in  j .  The 
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precision  in  the  determination  of  Mg  appears  a  little  better  than 
1  percent.  Therefore,  taking  into  account  the  variations  in  the 
individual  probe  readings  and  the  error  in  the  Mg  determination, 
the  accuracy  of  the  measurements  is  ±15  percent. 

As  shown  in  figure  6,  additional  measurements  were  made 
on  the  0.06-in.  electrode,  0.25-in.  diameter  probes  over  the  Mach 
number  range  7.2  to  11.0  at  probe  potentials  of  -21.6,  -10.8  and 
-5.4  v. 


All  measurements  with  the  other  two  sets  of  probes  were 
made  at  =  -  1008  v.  As  indicated  in  figure  7,  the  repeatability 
appears  to  be  about  the  same  as  previously  noted  in  figure  5. 

The  current  measurements  for  the  0,.02-in.  and  0.06-in. 


M 

s 

j 

0.02-in. 

electrode 

(|iamp) 

J 

0.06-in. 

electrode 

(|j.amp) 

j (0 .06  in.) 
j (0 .02  in.) 

8,6 

28.9 

245 

8.48 

9.0 

42.5 

355 

8.35 

9.4 

56 .1 

475 

8  46 

9.8 

69.7 

620 

8.89 

electrode  probes  with  diameter  of  0.25  in.  at  various  Mg  are  com¬ 
pared  in  the  above  table.  Assuming  no  fringe  or  sheath  effects  and 
no  variation  in  the  current  tangential  to  the  probe  surface,  the  j 
ratio  for  the  two  sets  of  probes  should  be  the  same  as  the  ratio  of 
their  electrode  areas,  9il.  As  shown  in  the  table,  this  ratio  is 
about  8.5:1;  the  difference  between  this  ratio  and  the  ratio  9:1  is 
within  the  accuracy  of  the  test  measurements.  A  reduced  ratio  could 
be  expected  to  occur  if  the  sheath  thickness  added  to^the  effective 
electrode  area,  for,  as  will  be  seen  in  section  4  (fig.  12a)  at  • 

(^p  =  -10.8  v,  6S  amounts  to  about  20  microns  over  the  Ms  range  8.6 
to  9.8.  Assuming  the  electrode  diameter  to  be  increased  by  26s, 
the  reduction  in  the  j  ratio  would  be  10  percent.  However,  since 
this  difference  is  also  within  the  experimental  accuracy,  fringe  or 
sheath  effects  may  not  be  accounted  for  in  this  manner.  Most  im¬ 
portantly,  however,  we  can  conclude  that  the  electrodes  are  located 
in  a  region  that  is  essentially  independent  of  variation  of  the  cur¬ 
rent  tangential  to  the  probe  surface.  Hence,  in  the  sheath  equations 
in  sections  3.7  and  4,  we  need  refer  only  to  the  one-dimensional 
Poisson  equation  for  the  electric  potential.  Although  the  measurements 
on  the  0.06-in.  electrode  probes  are  equivalent  to  those  of  the  0.02-in. 
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electrode,  0.25-in.  diameter  probes,  the  noise-to-signal  ratio  on 
the  latter  probes  was  significant  at  low  Mg.  Hence,  the  measure¬ 
ments  on  the  smaller  electrode,  0.25-in.  diameter  probes  were  dis¬ 
continued  and  future  reference  will  be  made  only  to  the  0.06-in. 
electrode  probes 

The  current  measurQments  shown  in  figure  7  for  the  0.22-in. 


M 

s 

j 

0.25-in. 
probe 
( |-iamp) 

j 

0 . 10-in. 
probe 
(pamp) 

J  (0.10  in.) 

J  (0.25  in.) 

8 

129 

48 

3.35 

9 

354 

154 

3.91 

10 

707 

305 

3.88 

11 

1288 

489 

3.42 

electrode,  0.10-in.  diameter  probes  are  compared  in  the  above  table 
with  those  shown  in  figure  6  for  the  0.06-in.  electrode,  0.25-in. 
diameter  probes  at  (|>  =  -10.8  v.  The  current  density,  J,  is  ob¬ 
tained  by  dividing  tne  current  by  the  electrode  area.  According  to 
the  table,  allowing  for  the  experimental  error,  the  measured  current 
densities  of  the  0.10-in.  probes  are  larger  by  a  factor  of  about  3.7. 
This  difference  is  due  to  the  fact  that  the  values  of  p  =  U/x  for  the 
0.10-in.  diameter  probes  are  larger  by  a  factor  of  2.5.  Since  the 
effect  of  (3  on  the  two  sets  of  measurements  will  involve  specific  cal¬ 
culations  of  boundary  layer  thickness,  gas  velocity,  ion-electron  re¬ 
combination,  etc,  the  discussion  of  the  effect  of  reducing  the  probe 
diameter  will  be  deferred  to  section  4. 1.1, (gee  pp.  30-32). 

In  the  solution  of  the  sheath  equations,  the  probe  poten¬ 
tial  must  be  referred  to  the  plasma  potential.  In  all  the  measure¬ 
ments,  the  probe  potentials  were  referred  to  the  shock  tube  wall, 
which  will  generally  be  at  a  potential  different  from  that  of  the 
plasma.  The  floating  potential  of  the  probes  was  found  to  be  +0.45 
±0.15  v  for  Mg  in  the  range  9.2  to  10.8.  In  gas  discharges,  the 
plasma  potential  is  approximately  0.2  v  larger  than  the  floating 
potential,  because  of  the  higher  random  velocity  of  the  electrons 
compared  with  the  ion  random  velocity  (ref  7).  In  the  boundary  layer, 
the  ion  random  motion  is  at  least  an  order  of  magnitude  larger  than 
the  uniform  gas  velocity.  Hence,  we  may  neglect  the  effect  of  the  gas 
motion  and  assume  the  plasma  potential  referred  to  the  shock  tube  wall 
as  +0.7  v. 


As  previously  noted,  a  fourth  set  of  probes  with  0.06-in. 
electrode,  0.25-in.  diameter  cylinder  was  constructed  having  an  outer 
brass  sleeve  concentric  with  the  cylinder  axis.  The  brass  and  glass 
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probe  measurements  were  compared  for  Mg  in  the  range  9.2  to  10.8 
at  (j>p  =  -21.4  v.  The  purpose  of  this  comparison  was  to  determine 
whether  a  static  charge  could  develop  on  the  glass  insulation  and 
thereby  influence  the  data.  Measurements  were  made  with  the  brass 
sleeves  both  insulated  from  and  grounded  to  the  shock  tube  through 
the  base  supports.  Since  within  the  experimental  error  of  ±15 
percent  these  measurements  appeared  identical,  only  the  glass  probe 
results  are  given. 

The  signal  forms  of  the  oscillograms  for  all  the  flat¬ 
headed  probes  appeared  approximately  the  same,  regardless  of  cylin¬ 
der  diameter.  The  upper  trace  in  figure  8  for  Ms  =  8.9  is  typical 
for  the  flat-headed  probes.  The  probe  current  is  measured  by  the 
vertical  scale  and  the  time  trace  is  left  to  right  at  20  psec  per 
division.  In  figure  8,  the  signal,  begins  at  the  point  (a)  and  is 
unsteady  for  about  20  psec  to  point  (c).  The  signal  then  remains 
at  an  approximately  constant  value  until  the  contact  surface  reaches 
the  probe,  whereupon  it  rapidly  falls  to  zero.  Probes  that  were 
not  clean  also  exhibited  an  interval  of  increasing  ion  current  sig¬ 
nal  prior  to  the  appearance  of  the  contact  surface,  the  total  rise 
being  as  much  as  20  percent  of  the  steady-state  signal.  The  dura¬ 
tion  of  the  steady-state  ion  current  was  always  at  least  30  psec 
and  increased  with  decreasing  Mach  number  to  about  120  psec  at  Mg  = 

7.2  . 

The  initial  unsteady  signal  lasted  for  about  16  psec  at 
Mg  =  11.0  and  lasted  as  long  as  40  psec  at  Ms  =  7.2.  For  a  given 
test,  individual  deviations  in  the  time  duration  of  the  initial 
unsteady  signal  for  the  four  probes  could  be  as  much  as  ±10  percent. 

As  suggested  by  the  strong  initial  pulse  shown  in  the  upper  trace 
of  figure  8,  the  unsteady  signal  would  seem  to  be  the  time  required 
for  the  flow  to  attain  a  steady  state,  rather  than  the  time  to  at¬ 
tain  a  steady-state  ionization  level.  The  strong  initial  ion  cur¬ 
rent  is  probably  due  to  the  detached  probe  shock  moving  out  into  the 
flow,  thereby  increasing  the  shock  strength  and  the  ionization  level. 

To  support  these  contentions,  the  conical  probes  described 
in  the  previous  section  were  constructed  to  determine  the  time  re¬ 
quired  to  attain  a  steady  signal  and  observe  the  signal  form.  As 
shown  in  the  lower  trace  of  figure  8,  no  large  initial  pulse  or  probe 
current  appears.  The  signal  begins  at  point  (a)  and  rises  linearly 
to  point  (b) .  Here,  the  absence  of  a  strong  initial  pulse  is  ex¬ 
plained  by  the  formation  of  an  attached  conical  shock  of  a  strength 
much  weaker  than  the  detached  normal  shock  of  the  flat-headed  probes. 
Moreover,  the  unsteady  signal  of  the  conical  probe  gradually  increases 
to  the  steady-state  value  at  point  (b)  in  about  6  psec  which,  at 
Mg  =  8.9,  is  approximately  four  times  longer  than  the  time  required 
for  the  shock  to  travel  the  0.142-in.  length  of  the  exposed  electrode 
(ref  8).  This  additional  time  may  be  required  to  fully  establish 
the  ion  current.  However,  the  unsteady  signal  for  the  next  14  psec 
to  point  (c)  is  not  understood.  The  unsteady  signals  shown  by  the 
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two  traces  of  figure  8  between  points  (a)  and  (c)  are  typical 
and  were  largely  repeatable. 

No  noticeable  change  in  the  test  data  appeared  because 
of  variation  in  T0  between  295  and  300°  K,  This  appears  reason¬ 
able  in  view  of  the  following  consideration.  At  pQ  =  1  cm  Hg  and 
Mg  between  9  and  10 ,  a  change  in  TQ  of  5°  K  amounts  to  a  change  in 
Te  of  80  to  95°  K  (ref  8) .  At  300°  K,  the  difference  in  Te  between 
Mg  =  9  and  10  amounts  to  900°'  K.  As  shown  in  figure  6  for  (j)p  =  -10.8  v, 
3  increases  from  0  355  ma  at  Ms  =  9  to  0.705  ma  at  =  10.  Since  the 
increase  in  n  is  due  primarily  to  the  increase  in  temperature  (ref  4) 
it  appears  that  the  5°  difference  in  T0  will  provide  about  a  10-per- 
cent  increase  in  j.  As  may  be  seen  from  figure  6  and  reference  8,  the 
variation  in  j  due  to  a  change  in  T0  from  295  to  300°  K  is  approximately 
the  same  at  other  Mg  and  (j)p0  Thus,  this  variation  in  j  is  within  the 
±15-percent  accuracy  of  the  measurements. 

Although  a  few  measurements  at  T0  <  295°  K  appeared  to  give 
lower  values  of  j ,  this  effect  was  not  sufficiently  investigated  to 
make  a  quantitative  statement. 

3  *  ION  FLOW  IN  THE  STAGNATION  REGION  OF  A  BLUNT  BODY 
3  - 1  Flow  Model 

It  is  convenient  to  describe  the  flow  over  a  probe  by  con¬ 
struction  of  a  model  in  which  distinct  flow  regimes  are  assumed.  It 
appears  possible  to  assume  the  following  three  distinct  flow  regimes0 

(1)  The  detached  probe  shock  and  outer  edge  of  the  thermal 
boundary  layer  form  the  boundaries  of  a  region  in  which  the  incident 
ion  number  density  is  increased  by  a  factor  of  more  than  100  and  at¬ 
tains  thermodynamic  equilibrium  upstream  of  the  edge  of  the  thermal 
bounda  ry  1 a  ye  r  u 

(2)  The  outer  edges  of  the  thermal  boundary  layer  and  the 
ion  sheath,  which  envelopes  the  probe  electrode  and  is  completely 
imbedded  within  the  boundary  layer,  form  the  boundaries  of  a  region 
in  which  the  usual  supersonic  stagnation-point  boundary  layer  equa¬ 
tions  for  a  nonionized  gas  mixture  apply  together  with  the  conserva¬ 
tion  equations  for  ions  and  electrons.  The  latter  equations  include 
terms  for  ambipolar  diffusion  and  formation-recombination  of  ion- 
electron  pairs.  Due  to  the  low  fractional  ionization,  the  boundary 
layer  equations  for  the  neutral  gas  are  not  altered  by  the  presence 
of  ions  or  by  their  additional  formation  or  recombination., 

(3)  Electrons  are  sharply  repelled  at  the  sheath  edge. 

The  neutral  gas  equations  together  with  the  ion  conservation  and  the 
one-dimensional  Poisson  equations  determine  the  ion  flow  within  the 
sheath. 


As  will  be  shown  by  numerical  calculations,  this  model 

appears  valid  for  Ms  >7,  the  value  for  the  lower  Ms  limit  depending 

on  d)  ,  D  and  p  , 

Tw*  o  r  o 
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The  calculations  of  sections  3.2,  3.3,  and  3.4  show  that 
the  ion-electron  transit  time  through  region  (1)  is  sufficiently 
long  to  attain  equilibrium  ionization  at  the  edge  of  the  boundary 
layer.  The  derivation  of  the  stagnation-point  boundary  layer 
equations  through  the  flow  regions  (2)  and  (3)  are  given  in  sections 
3.5,  3.6,  and  3.7.  Also  in  3,7,  the  electron  penetration  of  the 
sheath  is  shown  to  be  negligible.  In  section  4,  the  numerical  cal¬ 
culations  of  the  boundary  layer  and  sheath  thicknesses  and  the  ex¬ 
cellent  agreement  between  the  calculated  and  measured  values  of  Jw 
indicate  the  validity  of  the  assumed  flow  model. 

3.2  Flow  Time  from  Detached  Probe  Shock  to  Outer  Edge  of 
-  Thermal  Boundary  Layer 

The  Mach  number  of  the  gas  flow  incident  to  the  probe  is 
supersonic  for  Mg  >  7.  Thus,  the  gas  temperature  and  other  fluid 
properties  change  appreciably  across  the  detached  probe  shock.  The 
calculated  gas  mixture  velocities  and  properties  in  the  different 
shock  tube  flow  regimes  were  determined  using  references  8  and  9, 
and  the  ionization  densities  were  determined  using  reference  4. 

The  temperature,  density,  and  pressure  behind  the  probe 
shock  differ  by  less  than  10  percent  from  their  stagnation  values 
for  p0  in  the  range  2.0  to  0.2  cm  Hg,  TQ  =  300°  K,  and  Mg  from  7 
to  10.  In  particular,  for  this  flow  regime,  the  difference  in  the 
ion-electron  densities  between  conditions  immediately  behind  the 
probe  shock  and  stagnation  amounts  to  about  35  percent  of  the  for¬ 
mer.  Consequently,  to  calculate  the  time  for  ionization  to  reach 
equilibrium,  the  region  between  the  detached  probe  shock  and  the 
boundary  layer  may  be  considered  one  of  constant  fluid  properties 
and  ionization. 

The  gas  flow  decelerates  approximately  linearly  in  the 
inviscid  region  behind  the  probe  shock,  and  the  velocity  at  the  edge 
of  the  boundary  layer  is  only  about  2  percent  of  the  velocity  immed¬ 
iately  behind  the  shock.  The  ion  velocity  is  the  sum  of  the  gas 
motion  and  the  relative  motion  due  to  ambipolar  diffusion.  Because 
of  ion  formation  between  the  probe  shock  and  the  boundary  layer, 
the  diffusion  velocity  is  in  a  direction  opposite  to  the  gas  motion. 
In  the  following  calculations,  the  diffusion  velocity  is  neglected. 

We  obtain  time  values  somewhat  smaller  than  those  actually  required 
for  the  flow  to  cross  the  region  bounded  by  the  outer  edge  of  the 
thermal  boundary  layer  and  the  probe  shock.  As  will  be  seen  in  sec¬ 
tion  3.4,  these  times  are  still  sufficient  to  attain  equilibrium 
ionization . 

For  the  linear  deceleration  of  the  gas  motion  in  the  invis¬ 
cid  region  behind  the  probe  shock,  the  approximate  time  for  the  ions 
to  traverse  the  distance  between  the  detached  shock  and  the  thermal 
boundary  layer  is  given  by 
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Table  II  summarizes  the  calculations  for  Mg  from  7  to 
10  at  pQ  =  1.0  and  0.2  cm.  Values  for  pQ  and  Te  were  obtained 
from  reference  8.  It  is  clear  from  these  results  that  for  sheath 
thicknesses  of  the  order  of  a  micron  or  larger,  Collisions  occur 
within  the  sheath  and  the  drift  velocity  and  probe  current  will 
be  governed  largely  by  ion  mobility  and  not  by  free-fall  diffusion. 
The  sheath  thicknesses  calculated  in  section  4  (fig  12a)  indicate 
that  this  is  indeed  the  case. 

3.4  Time  Required  to  Attain  Equilibrium  Ionization 


As  noted  in  section  3.2,  the  flow  in  the  region  between 
the  detached  probe  ;  shock  and  thermal  boundary  layer  has  approxi¬ 
mately  constant  values  of  temperature,  density,  pressure,  and  equi¬ 
librium  ionization  (if  attained)  equal  to  those  occuring  at  stag¬ 
nation.  The  time  required  for  the  gas  to  reach  equilibrium  ioniza¬ 
tion  for  stagnation  conditions  is  given  below  for  Mg  from  7  to  10  at 
pQ  =s  1.0  cm  and  Mg  =  9  and  10  at  pQ  =  0.2  cm. 

The  production  of  ions  and  electrons  at  the  fractional 
ionization  levels  encountered  here,  less  than  10"4,  is  governed  by 
the  reaction 


N  +  0  2  NO+  +  e  (1.1) 

At  fractional  ionization  concentrations  less  than  10“4,  the  electron 
and  NO  concentrations  are  approximately  equal. 

The  formation  of  equilibrium  values  of  dissociated  N  and  0 
takes  place  within  several  hundred  collisions  (ref  5,  12,  13).  The 
mean  free  paths  of  the  ions  and  neutrals  behind  the  probe  shock  are 
approximately  equal  to  the  values  of  L  given  in  table  II.  Since  the 
ion  and  neutral  random  velocities  are  about  2  x  (10e)  micron/sec,  the 
required  time  for  equilibrium  dissociation  is  less  than  0.1  jisec. 

This  time  is  negligible  compared  with  the  values  given  in  table  I  for 
the  time  required  for  the  ions  to  traverse  the  region  bounded  by  the 
probe  shock  and  thermal  boundary  layer.  Hence,  the  time  required  for 
dissociation  will  be  neglected  in  the  following  calculations. 

The  change  in  the  ion  or  electron  number  density  is  given 
by 

W/m*  =  dn/dt  =  Kf  <N>  <0>  -  Krn8  (3.3) 


where  the  formation  and  recombination  rate  constants  are  given  by 
equations  (62-b)  and  (64-b)  of  reference  5  as 


K{  =  5(  10”1 1  )  T"^  exp(-32 , 500/T) 


and 


K  =  3(  10“3  )  T"1'6 
r 


cm3 /ion-sec 


(3.4) 
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The  equilibrium  ion  number  density  increases  across  the 
probe  shock  by  a  factor  exceeding  100.  Consequently,  the  reduction 
in  dn/dt  by  recombination  is  seen  by  equation  (3.3)  to  depend  on  na , 
and  therefore,  recombination  becomes  significant  only  as  n  ng . 

Hence,  we  may  neglect  recombination  and  assume  that  the  time  required 
to  attain  equilibrium  ionization  is  governed  by  the  equation 

tf  =  ne/(Kf  <N>  <D»e  (3.5) 

The  calculations  for  t^  are  summarized  in  table  III,  where  the  values 
for  the  equilibrium  densities  of  <N>,  <P>  and  n  were  determined  from 
reference  4. 

A  comparison  of  the  times  calculated  in  tables  I  and  III 
shows  that  equilibrium  ionization  is  attained  when  Mg  and  pQ  are 
sufficiently  large.  Accordingly,  for  the  0.25-in.  probes,  it  ap¬ 
pears  that  equilibrium  ionization  occurs  ahead  of  the  thermal  boundary 
layer  at  Ms  >7.4  when  p^  >1.0  cm  and  at  Mg  >  8  3  when  p^  >0.2  cm. 

As  noted  in  section  3.2,  taking  into  account  the  decreased 
shock  detachment  distance  and  boundary  layer  thickness,  the  ion  transit 
time  from  the  shock  to  the  thermal  boundary  layer  for  the  0.10-in. 
probes  will  be  0.35  times  as  large  as  that  for  the  0.25-in.  probes. 
Hence,  for  the  0.10-in.  probes,  equilibrium  ionization  occurs  ahead  of 
the  thermal  boundary  layer  at  Ms  >7.9  when  pQ  >1.0  cm  and  at  Mg  >8,8 
when  pQ  >0.2  cm. 

In  section  4,  equilibrium  ionization  will  be  assumed  at  the 
edge  of  the  thermal  boundary  layer  for  Mg  >  7.43  at  p  =  1  cm  for 
both  the  0.25-  and  0.10-in.  probes.  For  the  latter  probes  at  low  Mg 
this  requires  an  increase  of  the  time  values  in  table  I "or  a  decrease 
of  the  time  values  in  table  III,  each  by  a  factor  of  2.  This  factor 
can  be  accounted  for  by  the  ion  diffusion  due  to  formation  between  the 
shock  and  the  boundary  layer,  which  reduces  the  ion  flow  in  the  direc¬ 
tion  of  the  gas  motion.  Also  this  factor  can  be  accounted  for  by 
assuming  a  somewhat  larger  value  for  K^,  for,  as  indicated  in  pages 
25  and  34-36  and  in  figures  13-17  of  reference  5,  the  rate  constants 
may  actually  be  three  times  larger  than  those  given  by  equation  (3.4). 

3.5  Axisymmetric  Laminar  Boundary  Layer  Flow  of  Dissociated  Air 
at  the  Stagnation  Point 


The  preceding  calculations  show  that  equilibrium  ionization 
is  attained  at  the  edge  of  the  thermal  boundary  layer  for  sufficiently 
large  DQ,  pQ  and  Ms .  The  fractional  ionization  increases  rapidly  with 
increasing  Mg ,  but  there  is  only  a  small  change  with  pQ  in  the  range 
2.0  to  0.2  cm.  In  section  4.1,  the  test  measurements  of  the  stagna¬ 
tion-point  ion  current  density  will  be  compared  with  the  numerical 
calculations  at  Ms  =  7.43,  9.88  and  10.43  for  pQ  =  1  cm.  Here,  the 
maximum  ion  number  density  fraction  is  0.98  (10“4)  and  occurs  at  the 
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edge  of  the  thermal  boundary  layer.  Ion-electron  recombination 
occurs  in  the  boundary  layer  prior  to  reaching  the  sheath,  and 
the  amount  of  recombination  increases  with  increasing  Mg.  Con¬ 
sequently,  as  will  be  seen  in  section  4.1,  due  to  recombination 
and  the  additional  velocity  due  to  ambipolar  diffusion,  the  frac¬ 
tional  ionization  at  the  edge  of  the  sheath  is  at  most  10“B  for 
Mg  <  10.43.  The  fractional  ionization  is  further  reduced  in  the 
sheath  by  the  ion  acceleration  caused  by  the  probe  potential.  Be¬ 
cause  of  the  low  ion-electron  concentration,  the  presence-  of  the 
charged  particles  is  assumed  to  have  no  effect  on  the  neutral  gas 
equations . 


In  particular,  the  ratio  of  the  energy  required  to  pro¬ 
duce  the  ionization  behind  the  probe  shock  to  the  air  mixture  total 
enthalpy  is  less  than  10“3  for  Mg  =  10.43  and  decreases  with  decreas¬ 
ing  M  .  Assuming  ion  formation  from  N2  and  02,  approximately  10  ev 
per  ion  are  required  (ref  5).  At  Mg  =  10.43,  the  corresponding 
energy  required  to  form  iig  s=  23.2  x  IO14  ion/cm3  is  3.72  x  10-3 
joule/cm? .  This  is  insignificant  in  comparison  with  the  air  mix¬ 
ture  total  enthalpy  of  8.89  joule/cm3,  calculated  with  the  help  of 
reference  4  and  table  IV.  Hence  the  initial  ionization  level  in  the 
shock  tube  upstream  of  the  probe  shock  has  no  effect  on  the  enthalpy 
behind  the  probe  shock. 

In  appendix  A,  Coulomb  forces  are  shown  to  have  a  negligi¬ 
ble  effect  on  the  diffusion  and  viscosity  of  the  dissociated  air 
mixture.  Consequently,  in  the  following  presentation,  the  stagnation- 
point  boundary  layer  flow  is  assumed  to  consist  of  a  binary  mixture  of 
neutral  gas  particles  and  ion-electron  pairs. 

Neglecting  the  presence  of  the  ion-electron  pairs,  a  pro¬ 
cedure  for  obtaining  solutions  of  the  dissociated  neutral  gas  motion 
in  the  boundary  layer  is  given  in  reference  14.  Since  these  solu¬ 
tions  will  be  used  to  obtain  solutions  of  ion-electron. flows,  some  of 
the  details  presented  in  the  Fay-Riddell  method  of  solution  will  be 
shown.  In  particular,  reference  will  be  made  to  their  Method  2  solu¬ 
tion  . 


The  laminar  boundary  layer  equations  for  axisymmetric  dis¬ 
sociated  air  flow  are  given  by* 


Continuity : 

(pur)  +  ( pvr)  =  0 

A  y 

(3.6) 

Momentum: 

PUUx  +  pvuy  =  -  Px  +  <MUy>y 

(3.7) 

Ene  rgy : 

puhx  +  pvhy  =  (kTy)y  +■  UPX  +  |j(uy)a  + 

rZV(hJ  -  :)CJy]y 

(3.8) 

♦  Subscripts  x  and  y  (or  T|  below)  are  used  to  denote  partial  differ¬ 
entiation  . 
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Species  Conservation:  puc .  +  pvc  .  =  (D.pc.  )  +  W,  (3.9) 

— - y  jx  y  jy  Jp  jy  y  j 

where  y  and  x  are  directions  normal  and  tangential  to  the  body  sur¬ 
face,  and  r  is  the  distance  of  the  probe  surface  from  the  axis  of 
symmetry.  For  the  flat-headed  circular  cylinder  probes  previously 
described,  r  =  x  in  the  stagnation  region.  The  enthalpy  of  the  mix¬ 
ture  is  given  by  h  =  £  c  (h.  -  h°)  .  The  terms  for  thermal  diffusion 
are  omitted  in  equations'3  (3^8)  attd  (3.9)  since,  for  stagnation-point 
supersonic  and  hypersonic  dissociated  flows,  the  motion  due  to  ther¬ 
mal  diffusion  has  been  found  to  be  negligible  compared  with  the  flow 
due  to  concentration  gradient  (ref  15). 

In  stagnation-point  flow,  because  of  symmetry,  all  depend¬ 
ent  variables  are  functions  of  y  only,  except  u,  which  is  proportional 
to  x  times  a  function  of  y.  This  includes  W j ,  since  the  formation  and 
recombination  of  atoms  are  functions  only  of  the  local  thermodynamic 
variables.  Consequently,  for  stagnation-point  flow,  the  exact  partial  dif¬ 
ferential  equations  may  be  transformed  to  exact  ordinary  differential 
equations . 


The  exact  ordinary  differential  equations  are  obtained  as 
functions  of  a  single  dimensionless  variable  7]  through  use  of  the 
Howarth-Mangler  transformation  equations 


T)  s 


U 

V2? 


frpdy 

o 


5  = 


rn  \i  Ur2dx 

”w  w 


(3.10) 


and  the  additional 
df/aT)  =  u/U 

l  =  oh/dwuw 


nondimensional  variables 


f 


fi 

I  (Sf/ftT).dTl 


o 

g  =  (h  +  ua/2)/he 


e  =  t/t 

e 


=  c ,/c . 

3  je 


where  the  subscript  e  refers  to  values  in  the  "free  stream"  at  the 
stagnation  point.  Here,  "free  stream"  denotes  the  fluid  at  the 
outer  edge  of  the  boundary  layer  that  has  passed  through  the  de¬ 
tached  shock  close  to  the  axis  of  symmetry. 

To  obtain  a  solution  of  the  continuity  equation  (3.6), 
we  note  that  U  =  £x  with  (3  =  constant.  Referring  to  equation  (3.10), 
we  have 
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I 


^  =  X»v/Pp  u  /2 


ww 

5  =  8n  u  x3 

jx  KMw  w 

a"11  Tly  =  pV2P/0wMw 

We  now  assume  pur  =  *  ,  pvr  =  -  tjf  and  ijr  =  */2f  f  ( T])  -  This  gives 

y  x  * 

u  =  £xf_  =  Uf„ 

t)  n 


(3.11) 


and 


Pv  =  -fV 2Ppw4w  -  WPpwnw/2)  xf^ 

Thus,  the  continuity  equation  (3.6)  is  satisfied. 

By  substituting  into  the  terms  of  equation  (3.7), 

puu  +  pvu  =  pP8x(f.®  -  2ff  ) 


nn' 


and 


-p  =  p  UU  =pp2x 
x  Me  x  He 

(pu  )  =  2p(32x(£f  ) 

y  y  7)7)  7) 


Consequently,  for  stagnation-point  flow,  the  momentum  equation  be¬ 
comes  (see  eq  A1  of  ref  14), 

«Vn +  ffnn +'  I  <pe/p  ‘  fn3)  = 0 

Before  transforming  the  energy  equation,  we  note  that 
Vh  =  [sc.(dh  /dT)]  VT  +  xXh.  -  h°)?c 


J 


and  we  set  c  =  7}c.(dh./dT)  =  TJc.c  . 

P  J  J  J  PJ 


Combining  equations  (3.8)  and  (3.9),  there  results  (eq  20  of  ref  14), 

c  (puT  +  pvT  )  =  (kT  )  +  up  +  p(u  )2  -  £W.  (h.-h?)  +  yc  .D.pc.  T 
PHx  M  y  yy  *x  K  y  j  j  j  ^  PJ  jM  jy  y 

Assuming  0  is  a  function  of  T)  only,  the  terms  in  the  above 
equation  can  be  written 

c  (pu9  +  pv9  )  =  -2c  Ppf 9 
P  y  x  M  y  P  M  7) 

<kVy  =  (kVN^  =  (VVa)^2pp 


n 
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up 


/T  =  *-(p  UU  )u/T  =  f-/T  -  0  at  the  stagnation  point 

:  e  Kexe  KeTfe 


4U 


=  ^  *  f__sTl,  3  =  ~r-  X, i-J  -O  at  the  stagnation  point 


™  V  Vw^w  Tit! 


Fc  D  pc .  9  =  — 0  )[£  .D  c  =  — 9  yc  .L.c 

pj  y  jy  y  pwpw  tt  pj  j  jn  o  n  pj  j  jTi 


where  L  /•?  =  D.p/u.  Combining  terms,  the  energy  ecuation  becomes 
(eq  A10Jof  refJ14) 

w . ( h . -h°) 

~  Q  '  ■  -  *Q  +  ^  9j£c  .L  cJm  -  >  -----  -■=  0 


(-  c  9  )  +  c  f0  +  ^  9JCC  .L  c 

a  P  11  T!  P  T]  a  TTPJ,  JT 


23pT 


(i  .13) 


In  the  same  way,  we  find  that  the  species  conservation 
equation  (3  9)  becomes  (eq  34  of  ref  14), 


VWt,  +  fcjT! +  Y2^  - 0 


(3  14.) 


Fay  and  Riddell  assume  dissociated  air  to  be  £  diatomic 
gas  composed  of  "air"  molecules  and  "air"  atoms  with  properties  properly 
averaged  between  oxygen  and  nitrogen.  Accordingly,  we  call  cA  the 
mass  fraction  of  atomic  oxygen  and  nitrogen.  Now,  only  one  equation 
of  the  form  (3.14)  is  needed.  Closed  form  approximations  are  given 
for  the  thermodynamic  properties  as  functions  of  9  and  c A.  We  then 
have  the  system  of  three  equations  (3.12),  (3.13),  and  (3.14)  in  9, 
cA  and  f . 


In  their  Method  2  solution.  Fay  and  Riddell  obtain  two 
alternate  forms  of  solution  for  c^: 

(a)  The  cA  is  in  thermodynamic  equilibrium 

(b)  The  c^  is  frozen  in  the  boundary  layer,  so  that  W.~0, 
For  the  equilibrium  boundary  layer,  the  W.  term  is  eliminated  between 
equations  (3.13)  and  (3.14)  and  the  resulting  equation  is  solved 
simultaneously  with  the  momentum  equation  (3.12)  and  an  ecuation  that 
gives  cA  as  a  function  of  9,  For  the  frozen  boundary  layer,  equa¬ 
tions  (3.12),  (3.13),  and  (3.14)  are  solved  with  =  0 . 

J 

The  required  boundary  conditions  are 


f (0)  =  f^(0)  =  0 

9(0)  =  0W 


Y”>  ■ 1 
©(*)  =  1 


c^(0)  —  0  for  catalytic  wall 


i-  <  a(®> 


c.  (0)  =  0  for  nonctita lytic  wall  J 
AT) 


"Ac 


(3.15) 


18 


For  the  equilibrium  boundary  layer,  numerical  solutions 
were  obtained  by  Fay  and  Riddell  on  a  digital  computer  using  an 
iterative  procedure  whereby  values  were  selected  for  £-*.(0)  and 
0j|(O)  until  the  required  conditions  at  "infinity"  were  “Tiiet,  i.e., 
f  '  1  and  9  -*  1 .  For  the  frozen  boundary  layer,  values  were  also 

sSlected  for  either  cA-(0)  (catalytic  wall)  or  cA(0)  (noncatalytic 
wall)  with  the  added  condition  cA(»)  =  cAe . 

A  summary  of  the  various  parameters  and  initial  conditions 
for  which  numerical  solutions  of  the  stagnation-point  dissociated 
boundary  layer  flow  have  been  obtained  by  Fay  and  Riddell  appears  in 
table  II  of  their  Avco  publication.  For  use  in  the  present  work, 
these  authors  made  available  three  sets  of  numerical  calculations, 
which  are  identified  in  their  table  II  as 


(a)  equilibrium  boundary  layer  with  0^  =-  0.0538 

(b)  frozen  boundary  layer  with  0W  =  0.0497 

(c)  frozen  boundary  layer  with  0W  =  0.0797 

The  last  two  calculations  assume  a  catalytic  wall  with  respect  to 
atom  recombination.  All  the  calculations  were  made  for  Tw  =  300°  K, 
L.  =  L  =  1.4,  and  <j  =  0.71  According  to  reference  8,  the  calcu¬ 
lations  in  (a),  (b) ,  and  (c)  for  the  above  0W  correspond  to  the  shock 
tube  Mach  numbers  9.88,  10.43,  and  7.43,  respectively. 


The  above  solutions  give  f,  0,  and  l  in  terms  of  7].  With 
the  help  of  reference  4,  the  compressibility  factor  can  be  found. 
Denoting  this  factor  by  Z,  we  can  evaluate  p  =  p  /RTZ.  From  equa¬ 
tion  (3.10),  there  results 


- _  (.1 

y  =  J  (1/0)dTi 


(3.16) 


Hence,  all  solutions  may  be  expressed  in  terms  of  the  distance  from 
the  wall,  y,  when  pis  known. 


By  table  I,  the  Mach  number  of  the  flow  incident  to  the 
probes  is  between  2.1  and  2.7.  For  this  Mach  number  range,  it  ap¬ 
pears  that  3  may  be  determined  to  within  ±15  percent  by  the  rela¬ 
tion  given  in  reference  16  (see  also  ref  17-20) . 


3  0.64  a  /D 

e  o 


(3.17) 


where  a  denotes  the  stagnation  value  of  the  sound  speed.  Using 
the  values  for  pe  and  pe  given  in  table  I ,  a0  was  evaluated  with 
the  help  of  table  IV  and  figure  6  of  reference  9. 


The  functions  f,  df/dJj, 
at  Ms  =  7.43,  9.88,  and  10.43  for  the 


Ji  and  P  are  shown  in  figure  9 
0.25-in.  diameter  probes.  By 
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equations  (3.16)  and  (3.17),  y  oc*/d  .  Accordingly,  the  above  func¬ 
tions  for  the  0. 10-in  probes  may  al?o  be  obtained  from  figure  9  by 
reducing  the  given  y  values  by  the  factor  *7 2.5. 

Shock  tube  experimental  verification  of  the  Fay-Riddell 
theory  with  respect  to  stagnation-point  heat  transfer  is  given  in 
reference  21  for  flow  Mach  numbers  between  2.3  and  5.  This  compares 
with  the  flow  Mach  numbers  2.1  to  2.7  for  the  test  measurements 
described  in  section  2.  As  shown  in  reference  14,  the  stagnation- 
point  heat  transfer  for  equilibrium  boundary  layer  flow  is  approxi¬ 
mately  equivalent  to  that  for  frozen  flow  with  identical  catalytic 
wall  (for  atom  recombination)  and  "free  stream”  conditions.  However, 
the  two  flows  result  in  different  distributions  for  f,  0,  and  For 
a  frozen  boundary  layer,  the  assumption  of  a  noncatalytic  wall  (for 
atom  recombination)  results  in  excessively  small  heat  transfer  (ref 
14,21)  and  an  examination  of  the  recombination  rates  Indicates  that  a 
noncatalytic  wall  can  be  realized  only  at  initial  shock  tube  pres¬ 
sures  below  1  mm  Hg  (ref  21). 

3.6  Ion-Electron  Boundary  Layer  Flow  outside  the  Sheath 


Outside  the  sheath,  due  to  ambipolar  diffusion,  the  de¬ 
scription  of  the  ion  flow  is  the  same  as  that  for  electrons.  Accord¬ 
ingly,  only  the  equations  for  ions  will  be  considered.  The  ion  con¬ 
servation  equation  has  the  same  form  as  the  species  conservation 
equation  (3.14)  and  may  be  written  as 

(  -  L  s  \  +  fs  +  — —  m  0 

\a  i  Tl/T)  7]  2Pciep 

where  s  =  c./c  .  Equation  (3.18)  applies  in  the  region  6  <y 
lie  s  — 

Two  boundary  conditions  are  required  for  the  solution 
equation  (3.18).  With  the  help  of  equation  (3.16)  we  can  refer 
flow  and  boundary  conditions  as  functions  of  y.  At  y  =  6_,  0  = 

At  any  point  in  the  boundary  layer,  there  is  an  infinite  number 
(s,  8  )  pairs  that  yields  s  *  1  at  y  *  6T.  The  y  component  of  the  ion 
current  density  is  J  =  yen,  where  n  »  n  s/P  and  outside  the  sheath  r. 
v=r-v  +  D  s  /s.  Accordingly,  since  outlide  the  sheath  J  is  a  function 
of  s  andaSy,  the  other  boundary  condition  is  determined  by  the  require¬ 
ment  of  continuous  s  and  J  across  the  sheath.  Hence,  we  require  the 
(s,  J)  pair  at  6S  that  satisfies  the  boundary  condition  at  infinity  in 
equation  (3.18)  and  is  compatible  with  the  remaining  boundary  conditions 
given  for  the  sheath  equations. 

To  obtain  numerical  solutions,  a  value  for  the  sheath  thick¬ 
ness  was  arbitrarily  selected.  Starting  at  this  value  for  y  =  equa¬ 
tion  (3.18)  was  solved  for  a  number  of  (s,  s  )  pairs,  obtaining  thereby 


(3.18) 


=  T 
of 

to  the 
s  =  1 . 
of 
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)|c 

a  continuous  curve  of  J  as  a  function  of  s  at  6  .  By  equations 
(3.16)  and  (3.18),  at  any  6g,  we  see  that  the  curve  sys(ss)  must 
have  a  negative  slope.  From  the  above  definition  of  J,  for  the 
curve  Jc.(si5)  to  have  a  negative  slope,  the  increase  in  the  term 
Das  resulting  from  a  reduction  in  s  must  be  larger  than  the  re¬ 
duction  in  the  term  vs.  For  each  value  of  6S  where  0  >  <J>  >  -  30  v, 
all  calculations  demonstrated  that  the  curve  Js(ss)  had  a  negative 
slope  over  the  required  range  of  sg .  As  will  be  seen  in  section  3.7, 
the  possible  solutions  for  the  sheath  equations  lie  along  a  single 
curve  for  which  <3Jg/dss  >  0  and,  if  a  solution  of  the  boundary  layer 
flow  exists,  this  curve  intersects  the  previous  curve  at  a  single 
point.  The  solutions  for  the  ion  flow  through  the  boundary  layer 
are  therefore  unique  and  single-valued  for  each  5g. 

Equation  (3.18)  was  solved  using  the  solutions  for  £,  T, 
f,  and  p  given  in  section  3.5,  constant  values  for  L^/o  -  2  and  3, 
and  the  values  for  (3  and  cie/m^  =  ne/p  in  table  IV.  W  is  given  by 
equations  (3.3)  and  (3.4). 

The  contribution  to  W  from  ion-electron  formation  is  de¬ 
pendent  on  the  boundary  layer  distribution  of  <OI>  and  <D>.  The 
equilibrium  dissociation  fraction  at  6T  varies  between  14  and  23 
percent  at  pQ  =  1cm  for  7.43  <  Mg  <  10.43.  However,  ion-electron 
formation  in  the  boundary  layer  does  not  appear  significant  except 
in  the  vicinity  of  6T,  since  the  slow  rate  of  ion-electron  recom¬ 
bination  and  the  sharp  decrease  of  Kf  with  decreasing  temperature 
allows  for  ion-electron  number  densities  greatly  in  excess  of  equi¬ 
librium.  As  noted  in  section  3.4,  atomic  nitrogen  and  oxygen  recom¬ 
bination  occurs  very  rapidly.  Hence,  equilibrium  values  for  <N> 
and  <©>  are  assumed  in  the  boundary  layer.  These  were  obtained  with 
the  help  of  reference  4  and  are  summarized  in  table  V  for  Mg  =  7.43, 

9.88,  and  10.43  at  p  =1  cm.** 

7  o 


To  reduce  the  number  of  calculations  we  set  s  =  0  at  various  values 
of  y  <  6g  and  solved  equation  (3.18),  determining  therefrom  J  as  a 
function  of  s  at  each  6g. 

*  * 

The  value  for  W  depends  on  the  rate  constants  given  by  equation 
(3.4)  and  the  values  for  <N>,  <0>  and  n.  At  y  =  n  =  ne  refers 
to  the  equilibrium  value.  Equation  (3.3)  does  not  result  in  a  pre¬ 
cise  value  of  W  =  0  at  y  =  6-j.  In  the  calculations  of  equation 
(3.18),  we  set  W  =  0  when  ion-electron  formation  exceeded  recombin¬ 
ations;  otherwise,  the  calculated  value  for  W  was  used.  Of  course, 
due  to  the  large  boundary  layer  temperature  gradient,  ion-electron 
formation  never  exceeded  recombination  except  possibly  fpr  very  short, 
distances  from  the  edge  of  the  boundary  layer.  The  small  error  in 
the  calculation  of  W  causes  small  variations  of  s  ana  J  in  the  viciu- 
ity  of  6t  which  diminish  rapidly  with  decreasing  y  and  were  alwavs 
negligible  for  y  =  6S- 
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Ion  Diffusion — The  evaluation  of  J  requires  the  determina¬ 
tion  of  the  ion  motion  relative  to  the  neutral  gas  due  to  ambipolar 
diffusion.  No  values  are  known  to  the  author  for  ambipolar  diffusion 
coefficients  of  N0+  ion-electron  pairs  in  air.  However,  because  of 
the  structural  similarity  between  N0+  and  N2,  the  ambipolar  diffusion 
within  the  boundary  layer  should  be  approximately  twice  the  value  for 
self-diffusion  of  N2 .  Computations  of  the  latter  are  given  in  refer¬ 
ence  22  for  the  temperature  range  1,000  to  15,000°  K.  For  comparison 
with  these  computations  and  in  order  to  obtain  values  at  lower  tem¬ 
peratures,  calculations  were  also  made  for  twice  the  binary  diffusion 
coefficient  of  NO  in  air,  using  the  Lennard-Jones  potential  given  in 
reference  23.  In  the  notation  of  reference  23,  it  was  assumed  that 


M-^  =  30  =  molecular  weight  of  NO 

Mg  =  28  =  molecular  weight  of  partially  dissociated  air 

T  =  absolute  temperature,  °K 

p  =  pressure,  atmospheres 

^12  =  ^al  +  a2^  ~  3,54  where  CTj  =  3.47  X  for  NO  and 

cjg  =  3.62  X  for  air 

€  /*  =N/(£1/ft)(€ g/*>  =  107.3°  K,  where  ^  119°  K  for 

NO  and  €0/*  =  97°  K  for  air. 


Substituting  these  values  into  equation (8. 2  -  44)  of  reference  23, 
the  value  for  twice  the  diffusion  coefficient  of  NO  in  air  is  given 
by 


D  =7.80  (10“b)  T1*5)/?  0  (3.19) 

8.  0 

where  ft  =  Q(^T/€12>  is  the  Lennard-Jones  potential  given  in  table  I-M 
of  reference  23  and  pe  is  the  value  for  the  pressure  in  the  boundary 
layer . 


Comparing  the  above  calculations  for  NO  in  air  with  the 
tabulated  values  for  self-diffusion  of  N2  given  in  reference  22,  the 
latter  values  are  15  percent  larger  at  6000°  K,  become  equal  at  2500°  K, 
and  become  5  percent  smaller  at  1000°  K.  Since  this  agreement  was  so 
close,  equation  (3.19)  was  used  to  evaluate  the  ambipolar  diffusion 
coefficient  for  the  entire  boundary  layer  (exclusive  of  the  sheath, 
where  electrons  are  absent) .  The  results  of  the  calculations  for  PeDe 
are  given  in  table  V  as  a  function  of  temperature,  Da  was  obtained 
as  a  function  of  T)  for  each  Mg ,  since  T  is  a  known  function  of  7]  and 
p  =  pe  =  constant. 

Moreover,  because  the  pressure  is  constant  in  the  boundary 
layer,  the  value  of  Li/a  =  Dap/p  is  a  function  of  temperature  only. 


However,  L  /<j  varies  only  from  2.68  to  2  35  in  the  temperature  range 
300  to  5000°  K.  Thus,  as  previously  stated,  we  may  assume  L^/a  = 
constant  in  the  solution  for  s  in  equation  (3.18). 

3.7  Ion  Flow  within  the  Sheath 


The  current  measured  by  the  probes  is  due  entirely  to  ions, 
since  the  mean  electron  energy  is  only  about  0.3  v  and  the  electrode 
potential  was  always  strongly  negative,  at  least  -6.1  v  with  respect 
to  the  plasma.  Due  to  the  negative  probe  bias,  the  electrode  is  en¬ 
veloped  within  an  N0+  ion  sheath.  On  account  of  the  high  ion  density 
and  the  low  electron  energy,  the  transition  region  between  the  sheath 
and  the  neutral  plasma  is  very  narrow.  As  will  be  shown  at  the  end 

of  this  section,  the  effect  of  electrons  on  the  ion  distribution  and 
electric  potential  within  the  sheath  is  very  small.  In  particular, 
ion-electron  recombination  within  the  sheath  is  negligible.  Accord¬ 
ingly,  the  conservation  equation  for  electrons  may  be  dropped. 

The  amount  of  ion  formation  occurring  within  the  sheath  is 
negligible  compared  with  the  ion  flow  entering  the  sheath.  This  is 
due  to  (a)  the  slow  rate  of  ion-electron  recombination  in  the  boundary 
layer,  resulting  in  ion  number  densities  entering  the  sheath  that  are 
considerably  larger  than  those  which  would  occur  in  an  equilibrium  flow 
and  (b)  the  low  ion-electron  formation  within  the  sheath,  because  of 
the  low  values  of  <N>,  <D>,  and  Kf  associated  with  the  sharp  reduc¬ 
tion  in  the  boundary  layer  temperature. 

Within  the  sheath,  the  electric  field  and  potential  change 
rapidly  from  zero  at  the  sheath  edge;  the  potential  decreases  mono- 
tonically  to  the  applied  potential  at  the  electrode,  all  potentials 
being  referred  to  the  plasma  potential.  However,  in  describing  the 
effect  bn  the  ion  velocity,  it  will  be  shown  in  section  4.1  that  the 
nonuniformity  of  the  field  may  be  neglected,  and  the  field  is  every¬ 
where  one  of  moderate  strength.  Accordingly,  the  ion  motion  will 
vary  linearly  with  the  electric  field.  We  now  proceed  to  derive  the 
relations  necessary  to  evaluate  the  electric  potential  and  the  ion 
current  density  in  the  sheath. 

Taking  into  account  the  diffusion  of  ions  due  to  the  pres¬ 
ence  of  the  elctric  field  (|)  in  the  sheath,  the  conservation  equa¬ 
tion  (3.9)  for  ions  becomes^ 

vs  /P  =  (D.s  /P  +  sK(|>  /P)  (3.20) 

y  i  y  Ty  y 

where  W  =  0  in  the  sheath,  s  =  c./c.  =  nP/n  ,  and  c.  =0  due  to 

symmetry  in  the  thermodynamic  properties  at  the  stagnation  point. 

In  equation  (3.20),  refers  to  diffusion  of  the  NO+  ion 
through  the  dissociated  air.  The  mobility  and  ion  diffusion  coeffi¬ 
cients  are  related  by 
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or 


(3.21) 


K  =  eD^/AT 
K  =  11,600  Di/T 

D  and  K  may  be  determined  as  functions  of  y  for  any  given  Mg 
since,  by  section  3.6  and  table  y,  Di  =  D  /2  and  D  and  T  may 
be  expressed  as  functions  of  y. 

Poisson’s  equation  for  the  electric  potential  is 

4yy  =  -  en/€o  0.22) 

As  in  section  3.6,  the  y  component  of  the  ion  current 
density  is  defined  by  the  equation 


J  =  ven  (3.23) 

where  now,  for  flow  within  the  sheath,  the  y  component  of  the 
average  drift  velocity  of  ions  is 

v  =  -  v  +  D,  s  /s  +  K(t>  (3.24) 

i  y  'y 

Equation  (3.24)  gives  the  relation  between  s  and  s  . 

J 

Combining  equations  (3.20),  (3.20)  and  (3.24),  we  obtain 

J  /en  +  s ( v/P)  =  0 

y©  y 

Substituting  the  relations  (3.11)  and  (3.16)  for  v  and  dy/dT),  there 
results 


J  =  20f  en  s/P  (3.25) 

y  7)  e 


Except  at  the  stagnation  point  where 


f  =  0,  J  >  0. 

n  9  y 


It  is  of  interest  to  show  that  the  same  relation  (3.25) 
will  be  obtained  if  we  start  from  the  ion  conservation  equation 
in  the  form 


(PiUr>x  +  <PiVr>y  =  0 


where, 
we  see 
to  the 
causes 


as  in  section  3.5,  r  =  x,  u  =  £xf  and  n  =0.  Accordingly, 
that  the  loss  in  J  on  approaching ''the  probe  electrode  is  due 
gas  mixture  flow  tangential  to  the  probe  surface,  which 
a  corresponding  current  outflow. 
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Finally,  Poisson's  equation  (3.22)  may  be  written 


1014 
+  8.85 


en  s 
e 


P 


0 


(3.26) 


suits 


The  boundary  conditions  applied  at  6e  are 

P 

(a)  continuous  s  and  J  _  * 

( 3.27 ) 

(b)  *  =  $y  =  o 

Eliminating  v  in  equations  (3.23)  and  (3.24),  there  re- 

s  —  (1/D. )  [jP/en  —  sK(j>  tVS*1  (3,^o) 

y  i  e  y 


The  three  equations  (3.25),  (3.26),  and  (3.28)  are  explicit  expres¬ 
sions  for  the  derivatives  (6  )  ,  J  ,  s  in  terms  of  6  ,  J.  s.  Hence, 
if  initial  values  are  given  ar one  endJof  the  interval,  the  deriva¬ 
tives  can  be  calculated  and  a  step-by-step  numerical  integration  pro¬ 
cedure  is  possible. 

The  sheath  equations  cannot  be  integrated  from  the  sheath 
edge  all  the  way  to  the  wall,  since  integration  in  this  direction  is 
unstable  and  leads  always  to  either  positively  infinite  or  negative 
values  of  s  prior  to  reaching  the  wall.  Hence,  it  is  required  to 
integrate  from  the  wall  to  the  sheath  erl^e,  and  a  matching  procedure 
is  necessary  to  satisfy  the  boundary  conditions  (3.27). 

The  boundary  conditions  at  the  wall  can  be  given  in  either 
of  two  forms.  For  given  Jw,  (k  and  sw,  assuming  for  instance  sw  =  0 
(catalytic  wall),  we  can  calculate  the  wall  values  of  ancl 

s  and  proceed  at  once  as  indicated  in  the  step-wise  numerical  solu¬ 
tion.  If  on  the  other  hand  Jw,  and  s  are  given,  for  instance 

s yyf  ss  0  (noncatalytic  wail),  equation  (3.28)  can  be  used  to  evaluate 
sw.  Then,  the  procedure  is  the  same  as  in  the  first  case. 


In  comparing  the  two  cases  (i.e.,  either  s  =  0  or  s  =  0) , 
it  was  observed  that  the  results  for  <{>  ,  J,  and  s  asWfunctions^of  y 
were  identical  everywhere  in  the  sheath,  except  in  the  vicinity  of 
the  wall  for  a  distance  of  about  1  percent  of  the  sheath  thickness. 
Apparently,  this  means  the  quantity  contained  in  the  brackets  on  the 
right  side  of  equation  (3.28)  became  small,  so  that  approximately 


s 


(3.29) 


Accordingly,  only  the  values  of  (j)  and  J  were  important.  Thus, 
the  integration  for  (j>  ,  J,  and  s  ¥rom  theWwall  to  an  assumed  value  for 
6S  can  be  uniquely  determined  by  the  initial  values  for  and  . 
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To  obtain  a  continuous  s  and  J  across  the  sheath,  the  fol¬ 
lowing  procedure  was  chosen.  Values  for  (s  ,  J  )  obtained  from  the 
solution  of  equation  (3.18)  were  used  as  boundary  conditions  in  the 
solution  of  the  sheath  equations,  where  now  the  latter  integration 
was  started  at  the  sheath  edge  and  in  the  direction  of  the  wall.  Al¬ 
though  integration  in  this  direction  is  unstable  and  can  only  be 
carried  out  for  ghort  distances  from  the  sheath  edge,  the  required 
(s  ,  J  )  pair  has  a  value  of  s  that  lies  on  the  s(y)  curve  separ¬ 
ating  Ihe  previously  mentioned3 two  groups  of  positively  infinite  and 
negative  s  solutions.  Hence,  this  procedure  was  very  convenient  for 
determining  the  (s  ,  J  )  pair  that  satisfied  the  boundary  conditions 
(3.27a).  S  S 


Since  the  change  in  J  through  the  sheath  is  small,  values 
of  J  slightly  smaller  than  Jg  can  be  used.  For  each  value  of  Jw, 
we  may  choose  tj>  and  integrate  the  equations  from  the  wall  to  the 
edge  of  the  she^?h.  By  properly  adjusting  (J>  and  J  ,  we  can  obtain 
Jg  and  <j>  g  =  0  at  y  =  6S  •  At  the  same  time,^This  procedure  automat¬ 
ically  insures  the  proper  value  of  s_ . 

s 

In  carrying  out  the  calculations,  =  0  was  assumed. 

Since  the  reduction  in  J  through  the  sheath  never  exceeded  5  percent 
of  Js,  except  for  one  or  two  adjustments  of  Jw,  the  matching  procer 
dure  was  reduced  almost  entirely  to  finding  the  (j^  for  which  <J>ys  =  0.* 

By  equation  (3.26),  (j)  varies  linearly  with  s  and  there¬ 
fore,  by  equation  (3.28)  an  increase  in  Js  results  in  an  increase  of 
ss .  Hence,  the  (sg,  Jg)  pairs  form  a  single  curve  for  which  dJs/dsgX). 
At  each  the  solutions  of  equation  (3.18)  resulted  in  a  single 

curve  Js(ss)  with  a  negative  slope  Consequently,  each  solution  of 
the  ion  flow  through  the  boundary  layer  is  single-valued  and  uniquely 
determined . 

Finally,  (j)^  can  be  determined  by  a  direct  integration  of 
(j»  ,  starting  with  (|>s  =  0.  Thus,  each  6g  gives  a  value  for  (|>w.  Hence, 
a^particular  value  of  (^w  can  always  be  found  by  adjusting  6  ,  since 
(j>w  decreases  with  decreasing  5g .  Accordingly,  the  curves  for  (}y,  J 
and  s  may  be  found  as  functions  of  y  for  a  given  (j>w.  Moreover,  using 
6g  as  a  parameter,  curves  for  Jw((j>w)  and  may  be  obtained  for 

°S>  <j>w  >  ~  24  v 

Electron  Penetration  of  Ion  Sheath — The  electron  penetra¬ 
tion  of  the  ion  sheath,  Ay*  may  be  determined  from  the  relation  of 
the  Debye  shielding  distance  (ref  24). 


*The  solutions  for  the  ion  conservation  equation  (3.18)  were  ob¬ 
tained  on  an  IBM  7090  computer  and  the  solutions  for  the  sheath 
equations  were  obtained  on  a  PACE  131 -R  analog  computer,  manufac¬ 
tured  by  Electronic  Associates  Inc.  The  results  of  each  computer 
were  partially  checked  by  making  several  complete  sets  of  calcu¬ 
lations  on  the  other  computer. 
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Ay  =:  (*T/47me8>i=  6  90  (T/n)£cm 


To  determine  whether  the  electron  penetration  is  significant,  we 
compare  it  with  the  total  sheath  thickness.  For  example,  at  Mg  =  7.43 
and  y  =  6^  =  36.5  microns,  figure  9a  gives  0  =  0.92.  Since  T0  =  3764, 

T  =  3460«  K.  As  will  be  shown  in  section  4.1  for  <jw  =  -11.5  v  (fig. 
14a),  the  corresponding  ng  =  2.34  (lOia)/cin3  .  Thus,  the  electron  pen¬ 
etration  is  about  7.6  percent  of  the  sheath  thickness.  Similarly,  for 
(}w  <  -2  v,  the  electron  penetration  amounts  to  a  small  percentage  of 
the  sheath  thickness  at  other  values  of  and  for  other  M_  >  7„43„ 

b  b 

The  electron  penetration  will  reduce  j(Jl  j  in  the  neighbor¬ 
hood  of  y  =  6g  and  will  thereby  also  modify  s.  '^However,  for  a  given 
(^w,  this  effect  is  essentially  equivalent  to  increasing  the  sheath 
thickness  by  an  amount  smaller  than  the  electron  penetration.  Thus, 
the  effect  of  the  electron  penetration  is  small  and  will  be  neglected 
in  the  determination  of  the  ion-electron  boundary  layer  flow.  In  par¬ 
ticular,  the  sharp  rejection  of  electrons  at  the  sheath  edge  results 
in  the  cessation  of  ion-electron  recombination  inside  the  sheath. 

4.  SOLUTIONS  OF  THE  ION-ELECTRON  BOUNDARY  LAYER  FLOW  AND  COMPARISON 
WITH  MEASUREMENTS  OF  THE  WALL  ION  CURRENT  DENSITY 

4.1  Exact  Numerical  Calculations 


4.1  1  Discussion  of  Results 


In  accordance  with  the  procedure  given  in  section 
3.7,  exact  numerical  solutions  of  the  equations  for  ion-electron  bound¬ 
ary  layer  flow  are  presented  for  the  condition  of  frozen  atom  flow 
(Wj  =  0  through  the  boundary  layer)  at  Mg  =  7.43  and  10.43  and  for 
equilibrium  dissociated  flow  at  Ms  =  9.88  The  calculations  were  car¬ 
ried  out  for  the  "free  stream"  and  wall  conditions  given  in  table  IV 
and  the  boundary  layer  values  for  Di  =  Da/2,  <N>  and  <0>  given  in 
table  V.  These  results  are  summarized  in  figures  10  through  14 
for  the  0.10-  and  0, 25-in.  probes  at  pQ  =  1  cm,  TQ  =  =  300°  K  and 

-!<!),<  -24  v. 


Calculations  are  also  included  to  give  the  effects 
of  uncertainty  in  the  values  for  Li/<j,  n0,  K^,  and  .  Although 
L^/ct  —  2  —  constant  was  generally  used  in  equation  (3.18),  as  men¬ 
tioned  in  section  3.6,  L^/j  varies  between  2  and  3  in  the  boundary 
layer.  Both  microwave  reflection  measurements  and  magnetic  probe 
measurements  of  dc  conductivity  to  determine  electron  density  in  a 
shock  tube  (ref  5,  6,  25)  have  indicated  equilibrium  values  for  n 
somewhat  larger  than  the  theoretical  values  calculated  from  refer¬ 
ence  4.  In  addition,  reference  5  indicates  that  the  rate  constants 
Kj  and  Kr  lie  between  the  values  given  by  equation  (3.4)  and  three 
times  these  values.  To  evaluate  the  effect  of  these  uncertainties, 
the  calculations  were  repeated  for  L../(j  =  3,  the  n0  values  given  in 
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table  IV  increased  by  a  factor  2,  and  the  Kf  and  Kr  given  by  equa¬ 
tion  (3. 4). each  increased  by  a  factor  3.  The  repeated  calculations 
are  presented  only  at  Mg  =  7.43  and  10.43,  since  the  results  at  Mg=s 
9,88  and  10.43  appeared  to  be  similar  because  of  the  near  quality 
in  1^. 


When  nes  instead  of  s  is  considered  the  variable 
in  equations  (3.25),  (3.26)  and  (3.28),  we  see  that  nes  as  a  func¬ 
tion  of  y  in  the  sheath  is  independent  of  variations  in  the  param¬ 
eters  L^/a,  n0,  Kf,  and  Kr.  According  to  equations  (3.23)  and  (3.24), 
the  sheath  solutions  for  J(y)  will  also  remain  unchanged.  However, 
by  equation  (3.18),  changes  will  occur  in  the  solutions  for  J(n0s)  and 
J(s)  outside  the  sheath.  Hence,  changes  in  the  sheath  calculations 
due  to  variation  of  the  above  parameters  are  brought  about  by  the 

alteration  of  the  boundary  values  for  (s  ,  J  ) . 

s '  s 

The  changes  in  (sg,  Jg)  can  be  determined  approx¬ 
imately  in  the  following  way.  The  values  of  s  aqd  J  as  functions 
of  y  can  be  obtained  by  integration  of  equation  (3.18),  beginning 
the  calculations  with  s  =  0  at  y  -  0.  As  noted  in  section  3.7,  the 
correct  boundary  conditions  (sg,  Jg)  are  determined  by  a  matching 
procedure  involving  the  sheath  equations.  Consequently,  the  solution 
of  equation  (3.18)  with  s  =  0  at  y  =  0  will  generally  not  give  the 
correct  boundary  values  (sg,  Jg)  for  the  solution  of  the  sheath  equa¬ 
tions.  However,  the  changes  in  (s,  J)  due  to  variations  in  L j/g,  n0, 
Kf,  and  Kr  will  be  approximately  the  same  as  the  changes  in  (sg,  Jg ) . 
In  addition,  we  draw  attention  to  the  fact  that  the  calculations 
showed  that  the  variations  in  the  curves  Jw(^w)  were  mostly  governed 
by  variations  in  Jg  and  much  less  affected  by  variations  of  a  similar 
amount  in  sg .  Here,  for  a  given  (|>w,  Jw  decreases  with  decreasing  Jg  . 

Equation  (3.18)  was  integrated  with  s  =  0  at  y  =  0; 
figures  10  a  and  b  show  the  resulting  curves  of  s  (dashed)  and  J 
(solid)  as  functions  of  y.  The  calculations  were  repeated  for  each 
of  the  following  conditions: 

(1)  no  variation  of  the  parameters  (V), 

(2)  L±/a  increased  from  2  to  3  (0), 

(3)  Kf  and  Kr  each  increased  by  the  factor  3( 

and 

(4)  n^  increased  by  the  factor  2  (O)  • 

As  indicated,  these  parameters  were  varied  one  at  a  time.  The  sym¬ 
bols  within  parentheses  identify  the  calculation  condition  for  the 
curves  given  in  figures  10  and  11.  In  figures  11  and  12,  the  open 


♦By  equation  (3.3),  in  order  to  have  W  w  0  at  y  =  6^,  Kf  was  increased 
by  the  factor  4  when  ne  was  increased  by  the  factor  2. 
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and  solid  symbols  refer  to  the  0.25-  and  0.10-in.  probes ,  re¬ 
spectively. 

It  is  important  to  keep  in  mind  that  the  ne  and 
rate  constant  effects  on  J  and  s  are  interrelated.  Thus,  if  re¬ 
combination  in  the  boundary  layer  is  small,  the  curves  s(y)  and 
J(y)  will  be  little  affected  by  variation  in  the  rate  constants 
and  the  value  of  n0  will  be  important.  In  the  other  extreme,  the 
two  parameters  have  interchanged  roles  when  the  boundary  layer 
recombination  is  large.  In  addition,  since  ion-electron  recom¬ 
bination  ceases  at  the  sheath  edge,  the  effect  of  a  variation  in 
the  rate  constants  increases  with  decreasing  sheath  thickness. 

According  to  figure  10a  for  Mg  =  7.43,  the  var¬ 
iation  in  n0  will  produce  the  largest  change  in  Jfl  and  thereby  the 
largest  change  in  the  «*w(^w)  curve.  Moreover,  as  shown  in  figure 
10a,  increasing  the  rate  constants  reduces  the  J  values,  the  per¬ 
centage  reduction  of  J  increasing  with  decreasing  6g ,  However, 
referring  to  figure  12a,  6S  >  20  microns  for  (^w  <  -  4v  and  B  = 

0.25  in.  Hence,  the  effect  of  varying  the  rate  constants  will  be 
significant  only  at  small  negative  values  of  <j>w.  Furthermore, 
according  to  figure  10a  we  may  expect  the  variation  in  L^/a  to 
have  only  a  moderately  small  effect  on  the  Jw(^w)  curve.  Thus,  we 
see  in  figure  11a  that  the  calculated  curves  Jw(<j>w)  are  in  fact 
changed  in  accordance  with  the  above  changes  shown  in  figure  10a. 
Apparently  then,  at  low  Mg,  the  ion-electron  recombination  in  the 
boundary  layer  is  small  and  the  uncertainty  in  the  calculations  is 

most  sensitive  to  the  value  of  n . 

© 

At  Mg  =  10.43,  we  see  by  figure  12  a  that  3  <  6S 
<  11  microns  when  -2  >  <J>  >  -24  v  and  B  -  0.25  in.  From  figure 
10b,  we  observe  that  IheWlargest  variat?on  in  J  is  due  to  varying 
the  rate  constants.  Consistent  with  the  previous  remarks,  there 
appears  to  be  no  effect  on  J  due  to  varying  in  figure  10b  for 
the  above  range  of  5g .  At  any  £  in  this  range,  it  is  also  seen 
that  the  s  values  for  the  curve  (V)  a?.'e  twice  those  of  the  curve 
O)  *  Since  n  =  nes/P,  this  means  that  the  n  values  of  the  two 
curves  are  identical.  Apparently,  since  recombination  is  a  function 
of  na,  at  Mg  ■=  10.43  we  have  a  case  of  large  boundary  layer  recombin¬ 
ation  so  that  the  above  two  sets  of  calculations  asymptotically  at¬ 
tain  equal  values  of  J  and  n.  To  a  lesser  extent,  the  larger  recom¬ 
bination  occurring  at  higher  Mg  is  a  result  of  encountering  larger 
variations  in  the  boundary  layer  temperature  and  thereby  larger  var¬ 
iations  in  Kjf,  <N>  and  <0>, 

For  Mg  =s  10.43,  as  for  Mg  =  7.43,  we  observe  in 
figure  10b  that  the  variation  of  L^/c  should  have  a  small  effect  on 
the  JW(<L>  curves.  Referring  to  figure  11c,  we  see  that  the  above 
changes  in  the  Jw(<jw)  curves  are  in  fact  as  indicated  by  figure  10b. 
Consequently,  at  high  Mg,  the  uncertainty  in  the  calculations  is  most 
sensitive  to  the  values  for  the  rate  constants. 
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By  equations  (3.11)  and  (3  16) ,  since  0  a  1/DQ, 
reducing  the  probe  size  from  0.25  to  0.10  in.  means  an  increase 
in  the  gas  mixture  velocity  and  a  reduction  in  the  boundary  layer 
thickness  each  by  a  factor  of  s/ 2.5.  Hence,  the  amount  of  recom¬ 
bination  will  decrease  with  decreasing  probe  size  because  of  the 
more  rapid  ion-electron  transit  through  the  region  6T  >  y  >  6g . 
However,  the  values  for  p0  =  pw,  T  and  Tw  will  remain“uncEanged. 
Therefore,  n(y),  s(y)  and  J(y)  increase  with  decreasing  DQ • 

The  calculations  showed  that  the  changes  in  the 
J  ((j>)  curves  for  the  0.10-in.  probe  caused  by  variations  in  the 
alaove  parameters  were  approximately  the  same  as  those  for  the 
0.25-in.  probe.  In  figure  11a  for  the  0.10-in  probe  at  Mg  =  7.43, 
we  show  only  the  calculations  assuming  (a)  no  variation  of  param¬ 
eters  and  (b)  a  simultaneous  increase  in  ne  by  a  factor  2  and  in¬ 
crease  in  the  rate  constants  by  a  factor  3.  In  figure  11c  for  the 
0.10-in.  probe  at  Ms  =  10.43,  we  show  only  the  calculations  assum¬ 
ing  (a)  no  variation  of  parameters  and  (b)  an  increase  in  the  rate 
constants  by  a  factor  3.  These  calculations  give  the  approximate 
limits  of  the  values  of  Jw((j>w)  due  to  uncertainty  in  the  values  of 
the  above  parameters . 

The  test  measurements  of  Jw  for  the  0.25-  and  0.10- 
in.  probes  are  also  given  in  figures  11a, b,  and  c.  In  each  figure, 
the  test  measurements  are  given  at  (j>w  =  -6.1,  -11.5  and  -22.3  v  for 
the  0.25-in.  probe  and  at  (j^  =  -11  5  v  for  the  0.10-in.  probe.  In 
each  figure,  a  dashed  curve  is  drawn  through  the  experimental  points 
for  the  0.25-in.  probe.  Evidently,  as  shown  for  Ms  =  7.43  and  10.43, 
the  test  measurements  fall  within  the  limits  of  the  uncertainty  in 
the  calculations  obtained  by  the  above  variation  of  parameters.  The 
calculation  giving  the  best  fit  to  the  experimental  data  at  Ms  =  7.43 
for  both  probe  sizes  corresponds  to  simultaneously  increasing  ne  by 
the  factor  2  and  increasing  the  rate  constants  by  the  factor  3.  At 
Ms  =  10.43,  the  best  fit  for  the  0.25-in.  probes  is  obtained  by  in¬ 
creasing  the  rate  constants  by  the  factor  3,  whereas  the  best  fit  for 
the  0.10-in.  probes  corresponds  to  the  case  where  the  parameters  are 
not  varied.  Within  the  experimental  error,  the  comparison  between  the 
test  measurements  and  the  calculations  at  Ms  --  9.88  is  the  same  as 
that  for  Ms  =  10.43.  Because  of  the  near  equality  in  Ms,  we  may  as¬ 
sume  that  the  discussion  for  Mg  =  10.43  also  applies  at  Ms  =  9.88. 

For  the  condition  where  the  parameters  were  not 
varied,  the  sheath  thicknesses  are  given  as  functions  of  in  figure 
12a  at  the  aforementioned  Ms  and  D0.  By  equation  (3  26),  for  a  given 
y  =  «s,  IW  will  increase  with  increasing  n  =  nes/P  Since  n  is  not 
in  equilibrium  in  the  boundary  layer,  n(y)  increases  with  increasing 
Ms.  We  have  already  seen  that  n(y)  is  greater  for  the  smaller  probes. 
Accordingly,  as  shown  in  figure  12a,  for  a  given  (j>w  we  should  expect 
6S  to  decrease  with  increasing  Ms  and  decreasing  D0.  Therefore,  as 
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shown  in  figure  12b  for  a  given  (j>w,  we  should  expect  (j>  to  increase 
with  increasing  Ms  and  decreasing  DQ. 

The  boundary  layer  thickness  for  the  0.25-in.  probes 
is  shown  in  figure  9  to  be  about  70  microns  and  independent  of  Mg 
Therefore,  the  boundary  layer  thickness  for  the  0  10-in.  probes  is 
smaller  by  a  factor  of  s[ 2 . 5  and  amounts  to  about  44  microns.  A  com¬ 
parison  of  figures  9  and  12a  shows  that  6S  <  6j>  for  (j>w  >  -24  v  at 
all  Mq  when  D,,.  =  0.25  in.  and  at  Mc  =  9.88  and  10.43  when  D  =  0.10 

SO  S  q 

in.  However,  at  Mg  =  7.43,  extends  fully  to  the  edge  of  the  bound¬ 
ary  layer  at  (j>w  =  -  24  v  when  D  =  0.10  in 

The  numerical  solutions  of  the  complete  boundary 
layer  flow  were  obtained  for  Ms  =  7.43,  9.88,  and  10.43  with  DQ  =  0,25 
and  0.10  in  at  (j>w  =  -11.5  v.  The  solutions  for  the  flow  outside  the 
sheath  are  given  in  figures  13  ae-f  and  the  solutions  for  the  flow' within 
the  sheath  are  given  in  figures  14a-f .  The  quantities  J/JT,  s,  n/ne, 
v/vrp  and  Vjj/vDg  are  given  in  figure.  13  as  functions  of  y  from  the 
boundary  layer  edge  to  the  sheath  edge  and  (j>/(j>w,  ^y/^vw'  n/ns>  s/ss 
and  Sy/SyS  are  given  in  figure  14  as  functions  of  y  from  the  sheath 
edge  to  the  probe  wall. 

As  can  be  seen  from  figure  13  and  as  was  previously 
noted,  the  reductions  in  J,  n/ne  and  s  as  functions  of  y  increase  with 
increasing  M_  and  increasing  D_  due  to  the  larger  ion-electron  recom- 
bination.  In  all  calculations,  we  can  observe  the  approximate  linear 
variation  of  v  with  y  and  the  rapid  increase  in  the  value  of  vD  in  the 
neighborhood  of  the  sheath  edge.  Accompanying  the  increase  in  vD,  a 
rapid  change  in  SJ/dy  near  the  sheath  edge  can  be  seen  and  here 
SJ /dy  0  . 


These  changes  in  vD  and  J  near  the  sheath  edge  are 
not  due  to  a  sudden  change  in  recombination  but  rather  to  the  pres¬ 
ence  of  the  electric  field  within  the  sheath.  The  electric  field  : 
causes  a  rapid  acceleration  of  the  ion  flow  within  the  sheath,  thereby 
changing  the  values  of  ss  and  Jg  .  In  turn,  the  changes  in  the  bound¬ 
ary  conditions  (3  27a)  requiring  continuous  s  and  J  across  the  sheath 
affects  the  entire  boundary  layer  flow.  For  example,  at  Mg  =  7.43, 

D0  =  0.25  in.  and  (j)w  =  -11.5  v,  figure  13a  gives  6g  =  36.5  microns, 
sg  =  0.229  and  Jg  =  2.28  raa/cm2 .  Integrating  equation  (3.18)  with 
s  =  0  at  y  =  0  gives  the  minimum  values  of  s(y)  that  can  be  f-  tained 
in  the  absence  of  an  electric  field.  This  was  done  in  figu.  10a  for 
Ms  =  7.43,  where  s  =  0.600  and  J  =  1.54  ma/cm3  at  y  =  36.5  microns 
For  this  example  ,  a  comparison  of  the  two  cases  shows  that  the  elec¬ 
tric  field  causes  a  reduction  in  sg  and  an  increase  in  Jg .  The  re¬ 
duction  in  sg  is  consistent  with  the  increase  in  Jg  for,  as  was  men¬ 
tioned  in  section  3  6,  the  solutions  of  equation  (3.18)  give  curves 
for  which  BJs/ds^  <  0. 

In  each  figure  14,  if  the  y  scale  is  normalized  by 
the  given  value  of  6S,  the  normalized  solutions  of  the  flow  within 
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the  sheath  $v/$yw,  $/$w,  n/nB,  t  /s  ®nd  s/s  as  functions  of  y/6f 
appear  to  bey independent  of  Mg  and  D0.  The  sheath  solutions  fbr  J 
are  not  given,  since  in  the  sheath  J  is  constant  to  vithin  5  per¬ 
cent  of  J  .  The  reason  for  the  similarity  in  the  normalized  solu¬ 
tions  may  be  seen  from  the  approximate  sheath  equations  given  in 
section  4.2 


and 


S  B 


JP/eneK(jiy 

1014  J 
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(4.1) 

(4t3>- 


Assuming  J  and  K  are  constant  in  the  sheath,  =  0  at  y  =  6S, 

and  <j>  =  (j)^  and  <j>  =  (J>w  at  y  =  0,  successive  integrations  of  equa¬ 

tion  ^4.3)  give 

♦/♦y*  =  (1  -  y/6g)1/to  and  <Ww  =  (1  -  y/6s?/3 

Now,  by  equation  (4.1),  s/sw  =  (1  -  y/fig)”1^3 .  In 
each  of  figures  14,  we  see  that  s w/ss  &  0.07.  However,  a  good  fit 
to  the  s/sg  curves  is 

s/s  =  0.14<l-y/6fl)‘V^  (for  0.1  <  y/6  <  0.98) 

s  s  » 

Furthermore,  since  n  =  n  s/P,  a  good  fit  to  the  n/n 
curves  is  8 

n/n  =  sPj}  P  ^  0.20  (l-y/6  )'***  (for  0.2  <  y/fi  <  0.96). 

S  S  8  S  S 


Near  the  wall,  the  factor  (Pg/P)  must  be  included. , 

Lastly,  by  differentiation  of  s/sg  and  assuming  s  /s 
55  0.012  y 

sy/sys  =  0.012  (l-y/Bs)-9®  (for  0.1  <  y/6s  <0.95) 

As  previously  noted,  we  find  in  figures  13  and  14  that 
n(y)  increases  with  increasing  Mg  and  decreasing  DQ.  Hence,  in 
agreement  with  the  previous  results  given  in  figures  12a  and  b, 
figures  14  show  that  6S  decreases  and  increases  with  increasing 
Mg  and  decreasing  DQ. 

4.1.2  New  Technique  for  Evaluation  of  Kf  and  Kr 

As  was  observed,  the  ion-electron  flow  behind  the 
detached  shock  is  sensitive  to  variations  in  a  number  of  parameters. 
Because  of  this  dependence,  it  may  be  possible  to  devise  new  tech¬ 
niques  for  improving  our  knowledge  of  some  of  these  quantities. 

The  following  consideration  illustrates  a  procedure  for  the  deter¬ 
mination  of  the  ion-electron  formation  and  recombination  rate  con¬ 
stants  . 

The  flow  model  given  in  section  3.1  assumes  <  6T - 
However,  a  similar  model  will  apply  when  Lp  >  6S  6j,  providing 
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the  ion-electron  and  neutral  gas  quantities  can  be  described  at  6g. 
This  information  can  be  obtained  from  a  knowledge  of  the  formation 
and  recombination  rate  constants  and  the  diffusion  coefficient  for 
ion-electron  pairs  in  a  manner  similar  to  that  for  which  n  and 
were  calculated  in  the  boundary  layer  in  section  3.6.  In  addition, 
v  can  be  found  hs  described  in  section  3.5. 

The  electric  field  effects  on  the  ion-electron 
flow  can  be  varied  by  changing  <)>w  without  altering  the  gasdynamic 
effects.  Here,  a  change  in  (j>w  results  in  a  change  in  6g  without 
altering  6^  or  altering  the  flow  of  the  neutral  gas.  In  addition, 
both  the  gas  mixture  and  the  ion-electron  flows  will  be  changed  by 
varying  Mg,  Dq  and  pQ.  The  latter  variations  will  alter  both  6rp 
and  6g  .  By  suitable  variations  of  Mg,  <j>w,  DQ  and  pQ,  the  ion-elec¬ 
tron  flow  can  be  altered  by  changing  only  the  location  of  the  sheath 
edge  over  a  specified  interval  behind  the  detached  shock  when  L  > 

v 


Ion-electron  recombination  ceases  for  y  <  6S  be¬ 
cause  of  the  absence  of  electrons  in  the  sheath.  Furthermore,  since 
the  equilibrium  ion  number  density  increases  by  a  factor  of  more 
than  100  across  the  probe  shock,  until  the  ionization  approaches  or 
exceeds  that  for  the  equilibrium,  recombination  will  be  small  in  com¬ 
parison  with  formation  in  the  region  Lp  >  y  >  6S «<  Also  recombination 
is  small  when  the  sheath  edge  is  located~sufTiciently  close  to  the 
detached  shock.  By  varying  and  thereby  6S,  the  resulting  changes 
in  J  can  be  made  highly  sensitive  to  changes  in  the  ion-electron  for¬ 
mation  rate„  Furthermore,  by. v  varying  Mg,  we  can  alter  the  temperature 
behind  the  detached  shock.  Here,  we  affect  the  ion-electron  formation 
since  Kf  varies  sharply  with  temperature.  The  equilibrium  constant, 
K0,  is  known  to  within  about  10  percent  (ref  13),  and  the  relation  be¬ 
tween  the  formation  and  recombination  rate  constants  is  given  by  = 
KeKr.  Consequently,  for  various  locations  of  the  sheath  edge  suffi* 
ciently  close  to  the  detached  shock,  we  can  determine  the  values  of 
the  ion-electron  formation  and  recombination  rate  constants  as  func¬ 
tions  of  temperature  which  will  bring  the  calculated  and  measured 
values  of  stagnation-point  ion  current  density  into  agreement  This 
procedure  may  result  in  a  more  accurate  determination  of  the  rate  con¬ 
stants,  which  are  presently  known  only  to  within  about  a  factor  3 
(ref  5). 


4.1.3  Validity  of  the  Assumption  vE  =  K(j>y 

As  noted  in  section  3.7,  the  calculations  are  based 
on  the  assumption  that  the  ion  velocity  due  to  the  electric  field  is 
given  by  vjj  =  K(j>y  This  relation  is  valid  only  in  the  event  certain 
restrictions  are  satisfied.  We  now  examine  the  numerical  results  to 
determine  the  validity  of  our  use  of  this  relation. 

The  linear  dependence  of  vE  on  (j)y  does  not  hold  for 
very  strong  electric  fields.  An  indication  of  moderate  field  strengths 
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is  that  v  «  c,  where  c  is  the  random  thermal  ion  velocity  (ref  26). 
Strong  electric  fields  produce  large  values  for  vg.  Within  the  ion 
sheath,  the  ratio  c/v  decreases  from  the  sheath  edge  to  the  wall. 

At  the  wail,  T  =  300°  K  and  therefore,  c  =  4.6  x  104  cm/sec.  At 
the  wall  v  Vg.  For  the  range  0  >  <}  >  -24  v  and  Mg  <10.43,  the 
largest  value  for  v£w  =  2.7  x  104  cm/sec,  where  (^yw  <  5.2  x  104  v/cm 
is  given  in  figure  12b  and  Kw  is  given  in  table  IV.  Consequently, 
for  the  above  (^w  and  Mg  range  and  with  the  possible  exception  of  the 
immediate  neighborhood  of  the  electrode,  c/v  »  1,  as  required. 

In  order  to  be  able  to  neglect  the  effect  of  the 
nonuniformity  of  the  electric  field  on  the  ion  mobility,  it  is  nec¬ 
essary  that  the  product  |(j)  L|  «  j(j>  J  ,  where  L  is  the  ion  mean  free 
path  (ref  26).  Except  in  neighborhood  of  the  sheath  edge,  equa¬ 
tion  (4.3)  gives  the  approximate  relation  within  the  sheath 

*  =  .  J 

"vy  8.85  K<(v 


The  values  for  L  and  K  at  the  wall  were  previously  calculated  and  are 
given  in  tables  II  and  IV,  respectively.  The  ratio  E/K  is  nearly  con¬ 
stant  through  the  sheath,  since  both  quantities  vary  essentially  as 
1/p.  Therefore,  through  the  sheath  the  quantity  (|>yy  E  varies  as  l/<j>y 
Based  on  the  J  and  (j>v  values  given  in  figures  11  and  12b,  at  the  probe 
wall  we  obtain 
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for  7,43  ^  Mg  <.  10.43  and  (j>  <  -2  v .  Except  for  the  region  bounding 
the  sheath  edge,  amounting  to~about  5  percent  of  the  sheath  thickness, 
the  results  of  the  exact  calculations  given  in  figures  14a-f  show  that 
(^v  >0**  w •  Hence,  except  in  the  neighborhood  of  the  sheath  edge, 


Hy,1! 
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Moreover,  since  both  the  region  of  the  possible  elec¬ 


tric.  field  nonuniformity  and  the  value  for  (j)  within  this  region  are 
small,  a  nonuniformity  of  the  field  appear in^he re  will  not  seriously 
affect  the  ion  mobility.  Therefore,  the  conditions  under  which  Vg  = 
K(j>y  is  assumed  are  fully  met  and  the  use  of  this  relation  in  the  pres¬ 
ent  calculations  is  fully  justified. 


4.2  Approximate  Closed  Form  Solutions 


In  section  4.1,  we  required  the  use  of  electronic  computer# 
to  obtain  exact  numerical  solutions  of  the  ion-electron  boundary  layer 
flow.  Here,  a  simple  procedure  is  given  to  determine  the  approximate 
ion-electron  recombination  in  the  boundary  layer  and  the  approximate 
values  for  J  and  n  at  the  sheath  edge.  In  addition,  the  exact  sheath 
equations  are  simplified  and  numerical  results  for  the  flow  within  the 
sheath  can  be  obtained  from  approximate  closed  form  solutions.  As 
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will  be  shown,  the  approximate  solutions  compare  very  well  with 
the  previous  exact  solutions  and  also  yield  values  for  Jw  that  are 
in  close  agreement  with  the  test  measurements. 

We  start  by  simplifying  the  exact  sheath  equations.  Firs  4., 
we  observe  that  vE  »  v  within  the  sheath,  except  in  the  neighborhood 
of  the  sheath  edge.  Then,  by  equation  (3.29),  the  approximate  rela¬ 
tion  for  s  within  the  sheath  is  given  by 


and 


s  =  JP/en0K<jly 

(4.1) 

equation  (4.1),  equations 

(3.25)  and  (3.26)  become 

Jy  =  2pfT) 

(4.2) 

.  101  4  j 

(-.3) 

fyy  “  "  8.85  IC^T 

From  equation  (4.2),  we  see  that  the  change  in  J  will  be  small  for  small 
sheath  thicknesses.  In  all  the  calculations  of  section  4.1,  the  re¬ 
duction  in  J  through  the  sheath  never  exceeded  5  percent  of  Jg.  Ac¬ 
cordingly,  we  may  neglect  the  variation  of  J  and  drop  equation  (4.2) 


In  order  to  integrate  equation  (4.3),  the  direction  of  y 
is  reversed  and  it  is  assumed  in  equations  (4.4)  and  (4.5)  that  y  =  0 
at  the  sheath  edge  and  increases  positively  to  y  =  6S  at  the  probe 
electrode.  Then  for  the  conditions  ((  =  d  =  0  at  y  =  0*  and  assuming 
J  and  K  are  constants,  successive  integrations  of  equations  (4.3) 
result  in  equations 


(4.4) 


(4.5) 


In  equations  (4.4)  and  (4.5),  at  the  probe  electrode 
y  =  6  ,  (1)  =  (b  ,  =  6  and  J  =  J 

J  S>  •  TW>  Ty  TyW  W 

Numerical  solutions  of  the  sheath  equations  can  be  obtained 
when  the  value  for  J  is  known,  since  we  can  solve  for  (L(y),  <j>(y)  and 
s(y)  by  using  equations  (4.4),  (4.5)  and  (4.1).  By  definition  and 
since  J  is  practically  constant  in  the  sheath 

J  =  J  =  v  en  (4.6) 

s  s  s 

♦At  y  =  0,  the  Condition  6=0  leads  to  infinite  <j>yy  in  equation 
(4.3).  Hence,  in  the  neighborhood  of  the  sheath  edge  where  vg  = 

K<L  is  small,  the  omission  of  the  terms  (v  +  vD)  results  in  increased 
values  for  6yyj  in  turn,  this  results  in  numerically  higher  values 
for  <L  and  (j>.  Elsewhere,  vE  »  (v  +  vD)  and  their  omission  is  neg¬ 
ligible.  Accordingly,  neglecting  v  and  vD  in  the  sheath  results  in 
slightly  larger  values  for  (j>w,  and  vw. 
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To  evaluate  v  ,  reference  must  be  made  to  the  exact  calculations 
given  in  section  4  1  Then,  taking  into  account  the  ion-electron 
recombination  in  the  region  6S  <  y  <  6T,  a  value  for  Jg  corres¬ 
ponding  to  a  given  (j>  can  be  obtaineB.  The  following  presenta¬ 
tion  gives  the  details  for  the  evaluation  of  Js. 

In  figures  13,  except  near  the  sheath  edge  where  the 
ion  flow  rapidly  accelerates,  it  can  be  observed  that  v  is  nearly 
constant  through  the  boundary  layer.  For  each  case,  it  can  be  as¬ 
sumed  that  this  value  of  v  ^  vg .  However,  if  possible,  it  would  be 
preferable  to  express  vs  in  terms  of  the  gas  mixture  velocity  and 
calculate  the  latter  from  equation  (3.11).  In  this  way,  an  approx¬ 
imate  method  for  evaluating  vg  could  be  established  without  having 
to  refer  to  the  exact  calculations.  Accordingly,  the  exact  calcu¬ 
lations  for  v  were  compared  with  the  approximate  values  for  Vij,, 
where  the  latter  were  calculated  from  equation  (4.8)  given  below. 

It  was  found  that  approximately 

vg  =  -2vt  (4.7) 


For  the  calculations  given  in  the  previous  section,  the  validity  of 

equation  (4.7)  appeared  to  be  independent  of  both  M_  and  D  . 

s  o 


To  evaluate  vT,  we  have  df/d7)  s  u/U  &  9 
Then  from  equation  (3.11) 

T| 


VT  =-(N/2Ppw(iw/pe)  J  PdTl 


P  from  figure  9. 


As  indicated  by  figure  1  of  reference  14,  and  in  reference  27,  P  ^  g, 
where  g  =  (h  +  u2/2)/h  and  h  is  the  enthalpy  of  the  mixture  By 
figure  5  of  reference  14  for  equilibrium  or  frozen  dissociated  flows 
and  as  assumed  in  reference  1  for  frozen  flow,  g  =  0.477]  is  a  valid 
approximation  for  over  90  percent  of  the  thermal  boundary  layer.  The 
limit  of  integration  is  then  T)  =  2  1  and 

VT  =  -  (4'8> 


Moreover,  as  will  be  used  later  in  this  section,  we  see  that  1/p 
a  T  a  7]  a*/ y7  approximately. 

Using  the  above  result  for  /pdTl,  the  relation  for  the 
boundary  layer  thickness  given  by  equation  (3.16)  becomes 

6_  =  H  /20 

T  0  Kw^w 
Me 

As  noted  in  section  3.6  ion  formation  in  the  boundary 
layer  is  small  compared  with  reduction  through  recombination 
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Then,  assuming  recombination  only  in  the  region  6a  <  y  <  6^,  f or 
Ky  =  constant,  the  integration  of  equation  (3.3)  yields 

ne/ns  =  1  +  Xrnet  (4.10) 

Here,  t  =  0  at  y  =  6^  and  n(t  =  0)  =  ne. 

The  approximate  ion  transit  time  through  the  region 
6S  <  y  <  is  given  by 

t  =  (6t  -  6s)/v4  (4.11) 

For  given  v"g  and  fiT,  the  simultaneous  solution  of  equa¬ 
tions  (4.5),  (4.6),  (4.10)  ahd  (4.11)  yields  approximate  Values 
for  J,  $s,  t  and  ns  for  a  prescribed  (j>w.  Then,  as  previously 
noted  for  0  <  y  <  6g,  equations  (4.4),  (4.5),  and  (4.1)  can  be 

solved  for  (j>y(y),  <j>(y)  and  s(y).  To  illustrate  the  procedure, 
solutions  for  6S,  Jw  and  will  be  given  for  the  0.25-in.  diam¬ 
eter  probes  at  p  =  1  cm,  Tq  =  Tw  =  300°  K,  (j>w  =  -11.5  v,  and  = 

7.43,  9.88  and  10.43.  In  addition,  the  approximate  sheath  solu¬ 
tions  are  given  for  (j>^(y),  (}(y)  and  s(y)  at  Mg  =  9.88. 

The  approximate  values  for  vg  and  6T,  shown  in  table  VI, 
were  calculated  using  equations  (4.7),  (4.8)  and  (4.9)  and  the 
"free  stream"  and  wall  conditions  given  in  table  IV.  Dividing  the 
vs  in  table  VI  by  2,  the  resulting  approximate  vT  agree  to  within 
10  percent  of  the  exact  values  obtained  in  section  4.1  and  given 
in  figures  13.*  As  can  be  seen  from  figures  9  at  the  6T  given  in 
table  VI,  the  exact  values  for  0  are  between  0.92  and  0.97  and  the 
exact  df/d7]  are  between  0.96  and  0.98.**  Thus,  equation  (4.9)  ade¬ 
quately  defines  the  edge  of  the  boundary  layer. 

As  will  be  seen  from  the  calculations,  6g  ^  0.25  6ij,.  As  shown 
earlier,  1/p  a  T  a  H  a  \Ty  approximately*.  Consequently,  Tig  0.5  Tip  an^ 
•ps%»  2p  .***  The  product  of  K  and  p  does  not  appear  to  vary  appreciably 

*  This  agreement  holds  also  for  the  0.10-in.  probes. 

**  Since  the  y  scales  in  figures  9  and  the  6«j*  in  table  VI  are  all  re¬ 
duced  in  value  by  the  factor  .5  in  going  from  the  0.25-  to  the  0.10- 
in.  probes,  the  above  values  for  0  and  df/3Tl  also  result  for  the  latter 
probes . 

***By  neglecting  changes  in  Z,  we  introduce  an  error  in  ps  not  exceed¬ 
ing  10  percent. 
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with  temperature  (ref  26).  Using  equation  (3.21)  and  the  values 
for  Da  -  2D^  given  in  table  V,  to  within  a  factor  of  1.6  it  may 
be  assumed  that  K  -  6.0/p  cm3/v-sec  in  the  sheath,  where  p  is  in 
kg/m3.  Since  1/p  CSsfy,  an  average  gas  density  within  the  sheath 
will  be  assumed  as  "p  =  (2pg  -f  pw)/3,  and  the  corresponding  value 
for  K  is  K  -  6  0/p. 


By  table  IV,  T0  is  between  3764  and  6036°K  for  7.43  g 
Mg  £  10„43.  Since  Ts  Te/2  and  reference  5  gives  Kr  =  3  x  10”3/ 
T1'6  ,  Kr  =  1CT8  cm3/ion-sec  was  assumed  for  the  entire  region  5g  < 

y  ^  6t- 


The  results  of  the  calculations  for  Jw  are  summarized 
in  table  VI  Within  the  uncertainty  in  the  values  for  ng  and  Kr, 
it.  was  previously  noted  in  section  4„1  that  the  exact  calculations 
gave  values  for  Jw  equal  to  the  experimental  measurements.  Assum¬ 
ing  the  measured  Jw  are  representative  of  the  exact  calculated 
values,  table  VI  shows  that  the  approximate  calculations  for  Jw  are 
smaller  and  within  30  percent  of  the  measured  values.  However, 
since  the  approximate  calculations  for  the  Jw  are  all  smaller,  and 
considering  the  approximate  nature  of  equation  (4.7),  the  above  com¬ 
parison  could  be  improved  by  increasing  the  coefficient  in  equation 
(4.7)  from  2  to  2.3.  In  this  way,  the  approximate  values  for  Jw 
would  be  within  .13  percent  of  the  test  measurements.  To  this  same 
end,  Kr  and  n  could  be ^modified  within  the  limits  of  uncertainty 
noted  in  section  4.1 

The  exact  solutions  of  6e  and  6  given  in  figures  12a  and 
b  for  (pw  =  -11  5  v  at  Mg  =•  7U43  differ  by  a?  much  as  40  percent  from 
the  approximate  values  given  in  table  VI.  However,  the  agreement  be¬ 
tween  the  exact  and  approximate  iralues  of  these  quantities  is  much 
closer  at  Ms  -  9.88  and  10  43  Moreover,  the  exact  calculations  for 
$s  are  larger  at  all  Mg . 

The  i cn- electron  recombination  in  the  region  6S  <  y  <6^, 
is  indicated  by  the  ratio  ne/nA „  As  shown  in  table  VI,  the  approxi¬ 
mate  values  of  ne/r.c  are  smaller  than  those  given  for  the  exact  cal¬ 
culations  „ 

Since  the  Jw  for  the  two  sets  of  calculations  are  nearly 
equal,  the  higher  ne/ng  for  the  exact  calculations  is  mainly  due  to 
the  previously  noted  rapid  acceleration  of  tae  ion  flow  near  the 
sheath  As  in  section  4.1,  here  we  also  find  ne/ns  and  therefore  the 
ion-electron  recombination  is  small  at  Mg  -  7.43  and  increases  sig¬ 
nificantly  with  increasing  M&  .  Assuming  the  values  for  Tg  and  ps  ob¬ 
tained  in  section  4.1,  the  corresponding  equilibrium  values  for  ng 
were  calculated  using  reference  4  As  shown  in  table  VI,  both  the 
exact  and  approximate  values  for  ng  are  larger  than  those  for  equi¬ 
librium,  and  these  differences  increase  sharply  with  increasing  Mg . 
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In  figure  15,  a  comparison  is  made  of  the  approximate 
(dashed  curves)  and  the  exact  (solid  curves)  solutions  for  <|)  , 

tand  s  through  the  sheath  for  Ms  =  9.88  with  D  =  0.25-in.  at 

=  -11.5  v3  Here,  the  curves  are  all  drawn  so  that  y  =  0  repre¬ 
sents  the  wall  In  the  approximate  calculations,  use  is  made  of 
equations  (4.4),  (4:5)  and  (4.1)  and  the  single  value  K  =K  =  0.980 
throughout  the  sheath  in  equations  (4.4)  and  (4.5).  For  the  cal¬ 
culation  of  s,  since  P  =  pe/p,  K  =  6.0/p,  and  J  =  31.0  ma/cma, 
equation  (4.1)  reduces  to  s  =  20.3/(j>yO  Except  near  the  sheath  edge 
where  (J>  — ♦  0  and  due  to  the  difference  between  the  exact  and  ap¬ 
proximate  values  of  5S,  the  two  sets  of  calculations  for  6,  (|>  and 
s  agree  to  within  a  factor  2. 

Finally,  we  can  compare  the  previous  results  of  the  exact 
calculations  for  all  Mg  and  (|>w  with  the  results  that  can  be  obtained 
from  the  approximate  sheath  equations  (4.4)  and  (4.5).  Since  J  ^ 
constant  in  the  sheath,  equations  (4  4)  and  (4.5)  give  (J)  Ct*fy  and  (j> 
a  y1,6  Approximately,  these  are  the  functional  dependencies  shown 
in  figures  14  when  allowance  is  made  for  the  variation  in  temperature 
and  therefore  K  through  the  sheath,  particularly  iifear  the  wall  where 
K  rapidly  decreases. 

Also,  as  may  be  seen  from  figure  6,  both  J  and  the  ratio 
(^Wa/J  increase  with  increasing  (j>w  and  therefore,  by  equation  (4.5), 

6g  increases  with  increasing  (|)w.  By  equation  (4.5).  doubling  <j>w 
provides  an  increase  in  6S  by  a  factor  less  than  2T3  s  1.59,  This 
result  is  in  agreement  with  the  exact  solutions  for  SgC^)  shown  in 
figure  12a. 

Furthermore,  combining  equations  (4  4)  and  (4.5),  we  see 
that  (j>y  ol  ( J^)1^3  .  According  to  figure  6,  0  <  dJw/d<j>w  <  1.  Hence, 
doubling  (j)w  will  result  in  an  increase  in  (j>yw  by  a  factor  between 
1.26  and  1.59.  This  result  is  in  agreement  with  the  exact  solutions 
for  4yw^4w^  given  in  figure  12b. . 

5.  SUMMARY  AND  SUGGESTIONS  FOR  FURTHER  STUDY 


A  model  of  the  ion-electron  flow  in  the  stagnation  region  be¬ 
tween  the  detached  shock  and  the  wall  of  a  Langmuir-type  probe  is 
proposed  in  section  3.1  for  ionized  air  at  ion  mass  fractions  less 
than  10”4.  Based  on  the  close  agreement  between  numerical  calcu¬ 
lations  and  experimental  shock  tube  measurements  of  ion  current  den¬ 
sity  at  the  stagnation  point  of  the  probe,  it  can  be  concluded  that 
the  proposed  model  gives  an  accurate  description  of  the  ion-elec¬ 
tron  flow  in  the  stagnation-point  boundary  layer  when  negative  poten¬ 
tials  are  applied  at  the  probe  stagnation  point.  Although  this  was 
not  done  in  the  present  work,  these  solutions  can  also  be  extended 
to  cover  the  entire  stagnation  region  behind  the  detached  shock, 
since  the  ion-electron  formation  and  flow  velocity  can  be  determined 
in  this  region  by  the  method  given  in  section  3.6. 
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The  flow  model  given  in  section  3,1  assumed  6S  <  and 
the  model  is  valid  for  restricted  ranges  of  Mg,  <j>w,  DQ  and  p  . 

The  present  study  is  well  within  these  limits  and  covers  the  ranges 
7  <  Ms  <11  and  -  2  >  (|>w  >  -  24  v  with  DQ  =  0„10  and  0.25  in.  at  pQ 
=  1  cm.  A  similar  flow  model  will  apply  when  L  >  6g  >  6-j>,  and  sol¬ 
utions  of  the  ion^electron  flow  can  again  be  obtained  in  the  entire 
stagnation  region  behind  the  detached  shock. 

In  the  shock  tube,  the  ion-electron  number  density  in¬ 
creases  across  the  detached  probe  shock  by  a  factor  exceeding  100, 
and  this  increase  does  not  depend  on  the  state  of  the  undisturbed  in¬ 
cident  ionized  flow.  Hence,  in  the  shock  tube,  the  ion-electron  flow 
ahead  of  the  probe  shock  cannot  be  determined.  However,  when  ion- 
electron  formation  and  recombination  due  to  the  presence  of  the  body 
does  not  greatly  exceed  the  number  density  in  the  undisturbed  flow, 
we  can  trace  the  ion-electron  flow  from  a  point  in  the  undisturbed 
incident  flow  to  the  stagnation  point  of  the  probe  An  example  of 
this  situation  can  arise  for  flight  at  high  altitudes. 

Although  the  present  work  refers  to  a  Langmuir-type  probe 
in  the  form  of  a  flat-headed  circular  cylinder  with  an  electrode  lo¬ 
cated  at  the  stagnation  point  of  the  incident  flow,  in  principle,  the 
probe  may  have  any  geometry  for  which  at  least  atione  poiftt’  (the 'probe 
electrode)  the  gas  mixture  and  ion-electron  flows  can  be  determined. 
Thus,  it  may  be  possible  to  determine  the  ion-electron  distributions 
in  the  stagnation  region  of  any  axially  symmetric  blunt  body  and  in 
the  boundary  layers  of  conical  bodies  and  flat  plates.  Accordingly, 
it  would  be  of  interest  to  repeat  the  above  calculations  and  shock 
tube  measurements  for  these  types  of  bodies,  thereby  contributing  to 
the  further  understanding  of  ion-electron  boundary  layer  flows. 

The  ion-electron  flow  behind  the  detached  shock  has  been 
shown  in  section  4,1  to  be  sensitive  to  a  number  of  physical  proper¬ 
ties,  such  as  Da,  Dj ,  K,  Kr  and  Kf .  Because  of  this  dependence,  it 
may  be  possible  to  extend  the  present  work  and  develop  procedures  to 
improve  the  accuracy  to  which  some  of  these  properties  are  known.  An 
example  of  such  a  procedure  is  given  in  section  4.1  for  the  determina¬ 
tion  of  the  ion-electron  formation  and  recombination  rate  constants. 

The  procedure  given  to  obtain  numerical  solutions  is  of 
special  interest*  since  the  one-dimensional  Poisson  equation  for  the 
electric  potential  is  joined  with  the  usual  gasdynamic  stagnation- 
point  boundary  layer  equations  for  the  dissociated  neutral  gas  mix¬ 
ture  ,  Since  the  ion  mass  fraction  is  less  than  10"'4,  the  latter 
equations  are  not  affected  by  the  presence  of  the  charged  particles 
or  by  their  formation  or  recombination.  Hence,  solutions  for  the 
electric  potential  and  the  ion-electron  number  and  current  densities 
can  be  readily  obtained. 

However,  at  higher  levels  of  ionization,  the  enthalpy  and 
the  flew  of  the  gas  mixture  will  be  affected  by  ion-electron  forma¬ 
tion  and  recombination.  In  addition,  for  sufficiently  large  ioniza¬ 
tion,  Coulomb  forces  will  increase  the  cross  section  of  ions  for 
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collisions  between  charged  particles  and  thereby  tend  to  modify 
the  values  of  E,  and  |i.  The  calculation  comparing  the  energy 
for  ion-electron  formation  with  the  gas  mixture  total  enthalpy 
are  given  on  page  15  and  the  calculations  for  the  effect  of 
Coulomb  forces  are  given  in  appendix  A.  From  these  calculations, 
it  appears  that  ion-electron  formation  and  recombination  effects 
on  the  gas  mixture  enthalpy  as  well  as  Coulomb  effects  on  the 
boundary  layer  diffusion  of  ions  and  electrons  will  become  signif¬ 
icant  at  ion  mass  fractions  exceeding  about  10"2  In  a  shock  tube, 
for  a  given  probe,  the  fractional  ionization  increases  primarily 
with  increasing  Ms . 

On  the  experimental  side,  the  conductivity  of  the  gas 
increases  Ttlth  increasing  Ms  and  the  gas  may  provide  an  additional 
conducting  path  between  the  probe  electrode  and  the  shock  tube  wall 
This  effect  may  result  in  erroneously  larger  values  of  the  probe 
measurements.  A  procedure  to  determine  when  the  gas  conductivity 
affects  the  measurements  was  given  on  page  6.  Hence,  both  the  cal¬ 
culation  and  experimental  procedures  will  require  modification  at 
large  ion  mass  fractions. 

Finally,  the  present  study  can  be  extended  to  other  gases 
where  the  ionization  process  leads  to  positive  ions  of  a  single 
type  and  electrons  only.  To  date,  no  attempt  has  been  made  to 
treat  flows  consisting  of  other  types  of  charged  particles. 
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APPENDIX  A.  THE  EFFECT  OF  COULOMB  FORCES  ON  L,  AND  \i 

In  the  preceding  determinations  for  L,  D^  and  p,  collisions 
between  neutral  particles  and  between  ions  and  neutral  particles 
were  based  on  particles  treated  as  hard  elastic  spheres.  For  suf¬ 
ficiently  large  ionization,  Coulomb  forces  will  increase  the  cross 
section  of  ions  for  collisions  between  charged  particles  and  thereby 
tend  to  modify  the  values  of  the  above  quantities.  The  effect  of 
Coulomb  forces  will  be  shown  to  increase  with  increasing  Ms  and  be 
greatest  at  the  edge  of  the  thermal  boundary  layer.  In  the  follow¬ 
ing  calculations,  reference  will  be  made  to  the  case  for  Mg  =  10.43 
for  the  corresponding  "free  stream"  conditions  at  y  =  6T  given  in 
table  IV.  As  will  be  shown,  the  effect  of  Coulomb  forces  on  the 
quantities  X,  D^  and  ji  may  be  neglected  for  Ms  <*  10.43. 

By  equation  (5-63a)  of  reference  28,  the  assumption  of  binary 
collisions  is  valid  when  the  dimensionless  parameter 


X  =  AT/e8(2ne)1^  »  1 

where  the  subscript  e  refers  to  conditions  at  the  edge  of  the  bound¬ 
ary  layer.  The  present  calculation  gives  X  =  21.6,  and  we  therefore 
assume  binary  collisions  in  the  following  calculations. 

Mean  Free  Path — By  page  149  of  reference  28,  assuming  only  col¬ 
lisions  between  the  charged  particles  and  Coulomb  forces,  the  ion 
mean  free  path  at  y  =  6T  is 

!7  =  (3/2  Zn  )_1 
e  —  e 

where,  by  equations  (5-20)  to  (5-22)  of  reference  28,  the  collision 
cross  section  is 


Z  = 


st 


=  a/rr  (e2/*r)3  in  a 


with 


At 

87m  e 
e 


For  our  case.  In  A  =  4.45,  Z  =  1.210  (10“13)  cm3  and  =5  1.260 
(10“4)  cm.  By  the  calculations  of  section  3.3,  which  are  summarized 
in  table  II,  the  ion  mean  free  path  with  the  neutral  particles  can 
be  expected  to  be  L  <  10”B  cm  at  Mg  =  10.43  for  pQ  =  1cm  and  to  be 
smaller  toward  the  wall  .  Hence,  is  governed  by  collisions  between 
the  ions  and  neutral  particles. 


due 


The  ratio  n:<Ng>  decreases  within  the  thermal  boundary  layer 
to  ion-electron  recombination  and  due  to  ion  acceleration  within 


45 


the  sheath.  Moreover,  this  ratio  decreases  with  decreasing  Ms, 
(table  IV).  Hence,  we  may  conclude  that  L  is  governed  by  colli¬ 
sions  between  ions  and  neutral  particles  for  the  entire  flow 
field  at  Mg  <  10.43. 

Diffusion  Coefficient — For  diffusion,  the  collision  cross  sec¬ 
tion  for  molecules  treated  as  hard  elastic  spheres  may  be  obtained 
from  equations  (5-19)  of  reference  28  as 


Z  = 


4  No¬ 


where  a  is  the  average  diameter  of  the  colliding  molecules.  As 
noted  in  section  3.6,  a ~  ^  and  therefore  Z  =  8.89  x  10~1B 

cm3.  The  "first  approximation  to  the  diffusion  coefficient"  is 
given  on  page  245  of  reference  29  as 

L  12J1  2  |_2MJ  <Ng^ 

whe  re 

H  =  m1m2/(m1  + 

m^  =:  mass  of  neutral  particle  =  4.65  x  10“23 ;  g 
m2  =  mass  of  N0+  ion  =  4. 99  x  10“2a  g 

For  the  "free  stream"  conditions  at  the  edge  of  the  boundary  layer 
given  in  table  IV  at  Ms  =  10.43,  there  results  [D12^l  =  2,08  cms/sec . 
This  compares  with  the  diffusion  coefficient  of  64?5/2 (15 .57)  =2.07 
cm2/sec,  obtained  from  table  V.  (In  table  V,  values  are  given  for 
ambipolar  diffusion  of  N0+  ion-electron  pairs  in  air,  which  were  as¬ 
sumed  to  be  equal  to  twice  the  binary  diffusion  coefficient  of  NO  in 
air,  see  page  22  ,  and  the  factor  15.57  =  pe) . 

The  effect  of  Coulomb  forces  on  D.^,  like  that  for  L,  depends  on 
the  product  <N  >Z.  Since  the  mean  free  path  between  charged  parti¬ 
cles  is  more  than  a  factor  10  greater  than  the  mean  free  path  between 
ion  and  neutral  particles,  we  see  that  the  effect  of  Coulomb  forces 
on  the  diffusion  coefficient  is  negligible. 

Viscosity  Coefficient — The  coefficient  of  viscosity  for  a  pure 
gas  is  given  by  equation  (10,21,1)  of  reference  29  as 

_  5  (flTm)2 

^  ~  2  (22) 

z 

(22) 

where  m  =  molecular  mass  and  Z  =  collision  cross  section.  Treat¬ 
ing  molecules  as  hard  elastic  spheres,  Z^22)  =  8  4H  a10a  =  1 . 778X10”1 4 

1Z 


4C 


cm3  for  the  neutral  particles.  Assuming  ions  subject  to  Coulomb 
forces  only,  equation  (5-22)  of  reference  28  gives  Z^22)  =  2  Z^11^ 

=  2.60  x  10“ia  cm3  for  the  ions.  Since  the  temperature  and  mass 
of  the  ions  and  neutrals  are  approximately  equal,  the  ratio  of  the 
viscosity  of  a  gas  consisting  of  ions  to  that  of  a  gas  consisting 
of  neutral  particles  is  about  1:150,  or  the  inverse  ratio  of  their 
collision  cross  sections. 

The  "first  approximation  to  the  coefficient  of  viscosity”  of  a 
gas  mixture  is  given  by  equation  (12.5,1)  of  reference  29  in  the 
following  form 


where,  for  the  notation  of  reference  29,  and  are  the 

above  "first  approximations  to  the  coefficients  of  viscosity"  of 
the  constituent  neutral  and  ion  gases,  respectively,  n  =  n  /n  , 

n21  =  nl  =  n2  =  ne'  and  where  A  and  E  are  Quantities 

depending  only  on  the  interaction  of  molecules  of  different  kinds. 
For  hard  elastic  spheres,  pages  172-3  of  reference  29  gives  A  =  0.4. 
By  page  230  of  reference  29,  the  quantity 

2 

E  =  -  (nx  +  n2 ) ( D11  +  m2^Di2*^l  =  2,53  X  10~3  &/cm-sec 

where  [Dj2^i  =  2*^8  cm3 /sec  is  the  previously  calculated  "first  ap¬ 
proximation  to  the  coefficient  of  diffusion”  of  the  ion  and  neutral 
particles.  Substituting  numerical  values  into  the  above  equation, 
there  results 

[\J1  =  [li^  =  8.76  x  10”4  g/cra-sec 

The  quantities  n  and  E  decrease  with  decreasing  Mfl,  so  that 
the  equality  remains  between  and  Moreover,  n21  de¬ 

creases  in  the  boundary  layer  due  to  ion-electron  recombination  and 
ion  acceleration  within  the  sheath  due  to  the  electric  field.  Also, 
E  decreases  in  the  boundary  layer  approximately  as  «/t.  Thus,  the 
effect  of  Coulomb  forces  on  the  viscosity  coefficient  is  negligible 
for  all  Mg  g  10.43, 
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TABLE  I 


FLOW  TIME  FROM  DETACHED  PROBE  SHOCK  TO  OUTER  EDGE  OF 
THERMAL  BOUNDARY  LAYER  FOR  0.25-IN. DIAMETER  PROBES 


M 

s 

po 

cm  Hg 

vi 

m/sec 

Mi=Vai 

cm 

V2 

m/sec 

Pe 

atm. 

kg/m3 

e 

10B/sec 

t 

10”asec 

7 

1.0 

2080 

2.14 

0.254 

4.78 

0.445 

1.14 

49  6 

17.8 

8 

1.0 

2400 

2  0  33 

0.227 

7.05 

0.540 

1.26 

46.8 

15.2 

9 

1.0 

2750 

2  o  52 

624 

10.10 

0.625 

1.41 

45.6 

13.1 

10 

1.0 

3120 

2  .66 

690 

13.80 

0.710 

1.56 

45.2 

11.4 

9 

0.2 

2750 

2  .60 

595 

2.16 

0.140 

1 .39 

95.7 

10.8 

10 

0*2 

3120 

2.74 

0.204 

661 

2.96 

0.157 

1.50 

96.6 

9.41 

TABLE  II 

ION  MEAN  FREE  PATH  WITHIN  BOUNDARY  LAYER 


M 

s 

Po 

cip  Hg 

1 

Pe 

atm 

T 

e 

OK 

<x> 

g  e 

101 B/cm3 

micron 

L 

w 

micron 

7 

1.0 

4.78 

3570 

0.985 

0.182 

0.0153 

8 

1.0 

7.p5 

4090 

1  n  27 

0.142 

b 10104 

9 

1.0 

10.10 

4770 

1  .56 

0.115 

0.00725 

10 

1.0 

13.80 

5670 

1.79 

0.100 

0.00530 

9 

0.2 

2.16 

4560 

0.348 

0.516 

0.0340 

10 

0.2 

2  96 

5420 

0.402 

0.446 

0.0246 
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TABLE  V 


EQUILIBRIUM  <N>  AND  <0>  AND  AMBI POLAR  DIFFUSION  COEFFICIENT 

M  =7.43 
s 

M  =  9.88 
s 

M  =  10.43 
s 

_ l _ - _ 

T 

P  D 
e  a 

<N> 

<0> 

<N> 

<0> 

<N> 

<0> 

JO3  °K 

cra2/sec 

1/cm3 

101 7/cm3 

1/cm3 

101 7 /cm3 

1/cm3 

^10l7/cm3 

63.9 

8.89  (10)17 

58.2 

55.4 

3.74  (10)17 

52.6 

4.04 

60.6 

5.0  ! 

1 

47.3 

1.43 

51.6 

1.54 

59.6 

4.8 

44.2' 

0.914 

50.6 

0.990 

57.4 

4.6 

41.2 

0.564 

47.5 

0.611 

53.6 

4.4 

38.2 

0.333 

43.0 

0.361 

48.3 

4.2 

35.4 

0.187 

37.4 

0.203 

41.5 

4.0 

32.7 

08.9  (10)14 

30.7 

0.107 

33.9 

3.8 

30.0 

48.9 

24.1 

52.9  (10)14 

26.4 

3.6 

27.4 

14.5  (10)14 

10.7 

22.2 

17.3 

24.0 

18.9 

3.4 

25.0 

5.99 

7.40 

9.20 

11.8 

9.95 

12.7 

3.2 

|  22.6 

2.20 

4.73 

3.37 

7.38 

3.06 

8.00 

3.0 

|  20.2 

0.706 

2.77 

1.08 

4.28 

1.17 

4.63 

2.8 

I  18.1 

0.191 

1  .47 

0.292 

2.25 

0.316 

.  2.44 

2  .6 

16.0 

42.1  (10)11 

0.690 

64.4  (10)11 

1.06 

69.6  (10)11 

1.14 

2.4 

14.0 

7.15 

0.283 

i  11.0 

j 

0.433 

11.9  ! 

1  0.468 

i 

2.2 

12.1 

0.874 

0.0971 

}  1.34 

j 

j  0.148 

1.45 

i  0.161 

| 

2.0 

10.3 

0.0698 

0.0266 

!  0.105 

0.0407! 

0.116 

|  0.0440 

1.8 

;  8.65 

i 

i 

i 

i 

t  ! 

1.6 

7.06 

; 

l 

1 

i 

l 

i 

1  i 

i  ; 

1.4 

5.67 

I 

1 

i 

i  ! 

« 

1 .2 

4.41 

i 

* 

t 

*  i 

i  i 

i 

1.0 

3.27 

i 

l 

i 

j 

! 

1 

0.8 

2.25 

i 

| 

j 

i  ! 

•  i 

0.6 

1.39 

f 

| 

! 

1 

0.4 

0.696 

; 

0.3 

0.420 

j 

> 
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TABLE  VI 


APPROXIMATE  SOLUTIONS  OF  THE  SHEATH  EQUATIONS 
FOR  THE  0.25-IN.  PROBES  AT  -  11.5  V 


M 

s 

K 

V 

s 

- 1 

5t 

j 

j 

t 

i 

i - ; 

*  i 

1  ♦ 

Tyw 

Approx 
(sec  4.2) 
J 

w 

Measured 

!  J 

w 

cm3 

1 

•  1 

!  /  ! 

v/sec 

cm/sec 

micron 

|  mi cron 

jlO-6  sec 

i  V/cra  S ma/cma 

ma/cm3 

7.43 

2.14 

2230  ! 

48.0 

21.9 

1.17 

7,960; 

2.73 

3.12 

9.88 

0.980 

2800 

!  45.4 

7.49 

1.35 

23,200 

31.0 

35.9 

10.43 

0.855 

2910 

45.4 

6.92 

1.32 

i 

25,000 

J _ 

34.1 

50.8 

i 

Exact  Solution 
(sec  4.1) 

Approx 
(sec  4.2); 

Equilibrium 

i 

M 

s 

n 

e 

Ts 

Ps 

n  /n 
e  s 

i 

n  /n 
e  s 

t 

n  /n 
e  s 

101  4/cm3 

°K 

kg/m3 

7.43 

0.0840 

3226 

0.592 

3.58 

1.098  ; 

6.14 

9  .88 

11.16 

3242 

1.424 

34.4 

16.10 

473 

10.43 

23.2 

2837 

1.907 

45.0 

31.7 

884-62 

Figure  1.  Installation  of  the  probes  in  the  shock  tube.  The  four 
probes  are  shown  mounted  in  the  2-in.  thick  plexiglas 
section.  A  second  plexiglas  section  shown  adjacent  to 
this  section  was  not  used  during  the  test  measurements. 

In  the  background  are  the  two  Tektronix  555  dual  beam 
oscilloscopes  and  the  miniboxes  housing  the  mercury 
cells  and  resistor. 
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882-62 

Figure  4.  Typical  flow  about  the  probes  in  air,  Ms=9  and  Po=0.5  cm  hg. 
Illumination  of  the  hot  gas  caused  the  film  exposure. 
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Shocfc  Tube  Mach  Number 


Figure  5* 


Probe  current  vs,  shock  tube  Mach  number  for  four  0,25  inch 
diameter  probes  with  0,060  inch  electrodes,  ■  -10,8  volts. 
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1014-62 

Figure  8.  Typical  oscilloscope  traces.  Upper  trace  is  typical  for 
the  flatheaded  probes  and  lower  trace  is  typical  for  the 
conical  probes.  Time  traces  left  to  right  at  20  psec 
per  division. 
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Figure  13c*  Variation  of  J,  s,  n,  v  and  outside  the  sheath  for 
«  -11*?  volts  on  the  0*25  inch  probes  at  Mg  ■  10.143* 

J_  «  986  ma/cm2,  n  »  2*32  (lO^/crn^,  Vn,  ■  2,U60  cm/sec. 


v_  »  8,190  cm/sec.  6  ®  7*77  microns. 
Ds  s 
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Figure  13d. 


Distance  From  Wall,  Microns 


Variation  of  J,  s,  n,  v  and  outside  the  sheath  for 
0W  =  -11.5  volts  on  the  0.10  inch  probes  at  =  7®U3* 

JT  «  U.38  ma/c m2,  ng  =  8oU0  (lO12)/^,  «  2,820  cm/sec 

VDs  "  cm/sec.  5  a  30.5  microns. 


Distance  From  Wall*  Microns 


Figure  13  e.  Variation  of  J,  s,  n,  v  and  vQ  outside  the  sheath  for 

**  -Ho 5  volts  on  the  0*10  inch  probes  at  M  **9,88. 

W  2  1^  7  ® 

J^**66U  m a/cm  ,  n  “1*116  (10  ?)/cm  ,  vT»3,li80  cm/sec, 

v_  »  12,600  cra/sec*  5  »  7*25  microns. 


Figure  13f. 
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Distance  From  Wall,  Microns 

Yuriatioi  of  J,  s,  n,  v  and  outside  the  sneath  for 
PT.  -  -11.5  volts  on  the  0.10  inch  probes  at  M  =10.U3. 
Jm  =  1.55  amp/ cm2,  n  ■  2.32  9  v  «3,880cm/sec 

v_.  «  11,900  cm/sec.  6  **  5*6l  microns. 
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Figure  31; b.  Variation  of  /,  fi  •  n,  s,  s  through  the  sheath  for  /  »  -11.5 

j  y  w 

volts  on  the  0.25  inch  probes  at  M  »  9*83.  «  20,700 

o  ® 

volt/ cm,  ng  «  3.26  (10  )/cnr,  s^  »  O.OlUl;,  sys  »  lU5/cm 
and  6  *»  9.85  microns. 


s/s  *Su/u 


Distance  From  Wall,  Microns 


Figure  2k  c.  Variation  of  n,  s,  sy  through  the  sheath 

for  *  -11.5  volts  on  the  0.25  inch  probes  at 

Mg  -  10*1;3«  $  ■  26,200  volt/cm,  sg  »  0.0086l> 

n  -  5#  18  (101  -  118/cm  and 

s  ys 

6  »  7*77  microns, 
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Arrrrmti  rnmiim  i  nxt  nr:  t  mmi  mrmimn  i  rrrt  n : ;  i  n 

0  2  k  6  a 


Distance  From  Wall,  Microns 

Figure  l|*e.  Variation  of  0,  n,  s,  sy  through  the  sheath 

for  »  -11  •  5  volts  on  the  0.10  inch  probes  at 

Mg  -  9.88.  0  -  28,000  volt/cm,  sg  «  0.0279, 

n  -  5.8 9  (10*®)/cm3,  s,ra  -  367/cm  and 
s  ys 

6  ■  7.25  microns, 
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